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I. INTRODUCTION 

Particle physics is at a crossroads. Behind us is the 
standard model (SM), the remarkably successful theory 
of all known elementary particles and their interactions. 
Ahead of us is an equally remarkable array of possibili- 
ties for new phenomena at the weak scale. Never before 
has an energy scale been so widely anticipated to yield 
profound insights, and never before have there been so 
many ideas about exactly what these insights could be. 
In this article, we review the current state of experimen- 
tal searches for supersymmetry, the most widely studied 
extension of the SM. 



A. Motivations for New Phenomena 

There are at present many reasons to expect new 
physics at the weak scale m wea k ~ 100 GeV — 1 TeV. 
Chief among these is the Higgs boson, an essential com- 
ponent of the SM that has yet to be discovered. At the 
same time, there are also strong motivations for new phe- 
nomena beyond the Higgs boson. These motivations in- 
clude naturalness, dark matter, and unification. 

1. Naturalness 

The physical mass of the SM Higgs boson is given by 
m 2 h = ml 2 + Ami ■ (1) 

where mP h 2 is the bare mass parameter present in the 
Lagrangian, and the quantum corrections are 

. 2 A 2 f A d 4 P A 2 a2 
Am ~ — - 2 \ -f ~ — ^A 2 , 2 

where A is a dimensionless gauge or Yukawa coupling, 
and A is the energy scale at which the SM is no longer a 
valid description of nature. Because Am 2 is proportional 
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to A 2 ( "quadratically divergent"), it is natural to expect 
the Higgs mass to be pulled up to within an order of mag- 
nitude of A by quantum corrections ([TJ [2] [31 H). Given 
that unitarity and precision constraints require rrih to be 
at the weak scale (J5)), this implies A < 1 TeV, and new 
physics should appear at the current energy frontier. Of 
course, the Higgs boson may not be a fundamental scalar, 
but in this case, too, its structure requires new physics 
at the weak scale (JS). For these reasons, naturalness is 
among the most robust motivations for new physics at an 
energy scale accessible to accelerator-based experiments. 

2. Dark Matter 

In the last decade, a wealth of cosmological observa- 
tions have constrained the energy densities of baryons, 
non-baryonic dark matter, and dark energy, in units of 
the critical density, to be flTJ) 

fl B = 0.0462 ±0.0015 
Udm = 0.233 ±0.013 (3) 
fl A = 0.721 ±0.015 . 

The non-baryonic dark matter must be stable or very 
long-lived and dominantly cold or warm. None of the 
particles of the SM satisfies these conditions, and so cos- 
mology requires new particles. 

Perhaps the simplest production mechanism for dark 
matter is thermal freeze out (0 [5] ED rTTjl . In this sce- 
nario, a new particle is initially in thermal contact with 
the SM, but as the Universe cools and expands, this par- 
ticle loses thermal contact and its energy density ap- 
proaches a constant. Under very general assumptions, 
this relic energy density satisfies 



where (av) is the dark matter's thermally-averaged anni- 
hilation cross section. It is a tantalizing fact that, when 
this cross section is typical of weak-scale particles, that is, 
av ~ Q ;2 / m wcak' w here m wea k ~ 100 GeV, fix is near the 
observed value of f^DM given in Eq. ^ . If thermal freeze 
out is the mechanism by which dark matter is produced 
in the early Universe, then, cosmological data therefore 
also point to the weak scale as the natural scale for new 
physics. 

3. Unification 

The SM is consistent with the observed properties of 
all known elementary particles. It also elegantly explains 
why some phenomena, such as proton decay and large 
flavor-changing neutral currents, are not observed. The 
latter fact is highly non-trivial, as evidenced by the in- 
tellectual contortions required of model builders who try 
to extend the SM. 



At the same time, the SM contains many free param- 
eters with values constrained by experiment, but not ex- 
plained. The number of free parameters may be reduced 
in unified theories, in which the symmetries of the SM are 
extended to larger symmetries. In particular, grand uni- 
fied theories, in which the SU(3) x SU(2) x U(l) gauge 
structure is extended to larger groups, are significantly 
motivated by the fact that the SM particle content fits 
perfectly into multiplets of SU(5) and larger groups (fj"2")l . 
potentially explaining the seemingly random assignment 
of quantum numbers, such as hypercharge. 

One straightforward implication of the simplest ideas 
of grand unification is that the gauge couplings of the SM 
must unify when extrapolated to higher scales through 
renormalization group evolution. The gauge couplings 
do not unify at any scale given the particle content of 
the SM, but they do unify at the value gu ~ 0.7 at 
Af GUT ~ 2 x 10 16 GeV if the SM is minimally extended 
by supersymmetry (SUSY) and the supersymmetric par- 
ticles are at the weak scale (E2 ED ED ED Q21>- This 
unification is highly non-trivial, not only because the cou- 
plings are now so precisely measured, but also because gjj 
is in the perturbative regime and Mgut is in the narrow 
range that is both high enough to suppress proton decay 
and low enough to avoid quantum gravitational effects. 
This unification is only logarithmically sensitive to the 
superpartner mass scale, and the degree of its success 
is somewhat model-dependent; see, e.g., the review by 
Raby in Ref. (fl"8|) . In conjunction with the previous two 
motivations, however, it provides still more evidence for 
new physics at the weak scale, and selects supersymme- 
try as a particularly motivated possibility. 



B. Experimental context 

There are two main areas where new phenomena could 
appear in particle physics. Deviations from SM predic- 
tions could show up in measurements performed with in- 
creasing precision. Examples are the anomalies observed 
in the forward-backward asymmetry in the production of 
bb pairs in e + e~ collisions at the Z peak (see Ref. (|19[1 . 
particularly section 7.3.5), or in the anomalous magnetic 
moment of the muon (see, e.g., the review by Hocker and 
Marciano in Ref. (HE))- Even if such anomalies receive ex- 
perimental confirmation at a sufficient significance level, 
their interpretation will however remain ambiguous, be- 
cause it will involve virtual contributions to the relevant 
amplitudes of yet undiscovered, therefore most likely very 
massive, particles. The alternative approach is to try to 
observe directly the production of these new particles, 
which is among the goals of the experiments at colliders 
operating at the highest possible energies. 

The Large Hadron Collider (LHC) at CERN will soon 
occupy the energy frontier. When it comes into oper- 
ation, pp collisions will take place at a center-of-mass 
energy of 10 TeV, and of 14 TeV later on. The instanta- 
neous luminosity will be raised first to 10 32 cm -2 s _1 and 
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progressively to 10 34 cm 2 s 1 . With the enormous data 
samples accumulated, the two general purpose experi- 
ments at the LHC, ATLAS flU] [2T) and CMS (gl |23J, 
will be in a position to explore in great detail the physics 
at the TeV scale. Since this is an entirely new domain, 
and since there are strong reasons to expect new phenom- 
ena at that scale, as advocated in the preceding section 
of this review, it may well be that ground breaking dis- 
coveries are made at the LHC, even after a short period 
of operation, once the detectors are properly aligned, cal- 
ibrated, and well understood. 

Until then, the most constraining results on searches 
for new phenomena at high energy have been or are still 
being obtained at LEP, HERA, and the Tevatron. Pro- 
viding a comprehensive account of such searches for su- 
persymmetry is the purpose of this review. 

The large e + e~ collider (LEP) at CERN operated from 
1989 to 2000. In a first phase (LEP1), the center-of-mass 
energy was set at or close to 91 GeV, the peak of the 
Z boson resonance. Four experiments, ALEPH (|2"H |2"5")1 . 
DELPHI (HD|57J|, L3 (gHD, and OPAL (HO) studied mil- 
lions of Z decays that allowed them to perform strin- 
gent precision tests of the SM. From the end of 1995 on, 
the energy was progressively increased (LEP2) to reach 
209 GeV in the center of mass during the last year of op- 
eration. Altogether, each of the experiments collected a 
total of - 1 fb" 1 of data, of which - 235 pb" 1 in 2000 at 
and above 204 GeV, the data set most relevant for new 
particle searches. 

At DESY, the HERA collider operation was termi- 
nated in June 2007. There, e ± p collisions were collected 
by two experiments, HI (30) and ZEUS ()3"T)l . at a center- 
of-mass energy of ~ 300 GeV. This was an asymmetric 
collider, with e ± and proton beam energies of 30 and 
820 GeV, respectively. An upgrade took place in 2001 
(HERA2), leading to higher luminosities than in the pre- 
vious phase (HERA1), and allowing operation with po- 
larized e beams. The data sets collected at HERA1 
and HERA2 with electron or positron beams altogether 
correspond to an integrated luminosity of ~ 0.5 fb _1 per 
experiment. 

Until the LHC comes into operation, the highest en- 
ergy collisions are provided by the Tevatron pp collider 
at Fermilab. During its first phase of operation (Run I), 
the center-of-mass energy was set to 1.8 TeV, and a data 
sample of ~ 110 pb _1 was collected by each of the two 
experiments, CDF <$2) and D0 ®. The highlight of 
that period was the discovery of the top quark in 1995. 
Major upgrades of the accelerator complex and of the 
two detectors took place for the second phase (Run II), 
which began in 2001. The center-of-mass energy was 
raised to 1.96 TeV, and the instantaneous luminosity was 
progressively increased to regularly approach or exceed 
3 x 10 32 cm~ 2 s- 1 in 2008. More than 5 fb" 1 of in- 
tegrated luminosity had been delivered by the Tevatron 
by the end of fiscal year (FY) 2008, and it is expected 
that another ~ 1.5 fb _1 of luminosity will be provided 
per additional year of operation. At the time of writing, 



running in FY 2009 is underway, running in FY 2010 is 
increasingly likely, and running in FY 2011 is kept as an 
option. 

All general purpose detectors at colliders share simi- 
lar features. A cylindrical "barrel" structure parallel to 
the beam axis surrounds the interaction region, and is 
closed by "end caps" perpendicular to the beam. The 
first elements encountered beyond the beam pipe are 
charged-particle detectors, with those closest to the in- 
teraction point benefiting from the highest spatial pre- 
cision. This tracking system is immersed in an axial 
magnetic field provided by a solenoidal magnet. Beyond 
the tracking system, electromagnetic calorimeters pro- 
vide electron and photon identification and energy mea- 
surement. These are followed by hadron calorimeters for 
the measurement of jet energies. Finally, track detectors 
are used to identify and measure the muons which have 
penetrated through the calorimeters and possibly addi- 
tional absorber material. 

Non-interacting particles, such as neutrinos, are de- 
tected by an apparent non-conservation of energy and 
momentum. In e + e~ annihilation, the missing energy 
and momentum can be directly inferred from a mea- 
surement of the final state particles, by comparison with 
the center-of-mass energy of the collision. In hadronic 
or ep collisions, the partons participating in the hard 
process carry only a fraction of the beam energy, and 
the beam remnants associated with the spectator par- 
tons largely escape undetected in the beam pipe. As a 
consequence, only conservation of the momentum in the 
direction transverse to the beams can be used, and the 
relevant quantity is the missing transverse energy Ifix, 
rather than the total missing energy. 

The mass reach in pp collisions at the Tevatron is ex- 
pected to be substantially larger than at LEP because of 
the higher center-of-mass energy. However, since the ini- 
tial partons participating in the hard process carry frac- 
tions X\ and X2 of the beam energy, the effective center- 
of-mass energy is only y/l = ^Jx\Xis. Because of the 
rapidly falling parton distribution functions (PDFs) as 
a function of those energy fractions, increasingly large 
integrated luminosities are needed to probe larger and 
larger \f% values. At HERA, furthermore, the center-of- 
mass energy in the eq collision cannot be fully used for 
new particle production, except in some very specific in- 
stances. This is in contrast to e + e~ or qq annihilation, 
and to gg fusion. As a consequence, the most constrain- 
ing results on new particle searches typically come from 
LEP and from the Tevatron. 

In the following, all limits quoted are given at a confi- 
dence level of 95%. 



II. SUPERSYMMETRIC MODELS AND PARTICLES 

Supersymmetry (SUSY) (|34l 135) l36|) is an extension of 
Poincare symmetry, which encompasses the known space- 
time symmetries of translations, rotations, and boosts. 



4 



As with the Poincare and internal symmetries, SUSY 
transforms particle states to other particle states. In con- 
trast to these other symmetries, however, SUSY relates 
states of different spin, transforming fcrmions into bosons 
and vice versa. None of the known particles can be super- 
symmetric partners of other known particles. As a result, 
SUSY predicts many new particle states. If SUSY were 
exact, these particles would be degenerate with known 
particles. Since this is experimentally excluded, if SUSY 
is a symmetry of nature, it must be broken. 

SUSY is the most studied extension of the SM be- 
cause it directly addresses several of the motivations for 
new physics discussed in Sec. [I] In supersymmetric theo- 
ries, the quadratically divergent loop contributions to the 
Higgs boson mass from SM particles are canceled by sim- 
ilar contributions from superpartners, ameliorating the 
gauge hierarchy problem ([571 [551 Supersymmetric 
theories also include excellent dark matter candidates, in 
the form of neutralinos (|40l 141)1 and gravitinos (|42l l43)l , 
that may naturally have the desired relic density. Finally, 
SUSY is strongly motivated by the hope for unifying 
forces, as it makes gauge coupling unification possible in 
simple grand unified theories (GUTs) CL3J QU E2 [16j EJ . 
It is important to note that all of these virtues are pre- 
served only if the superpartner mass scale is around the 
weak scale. The existence of SUSY in nature, although 
not necessarily at the weak scale, is also motivated by 
string theory and the beautiful mathematical properties 
of SUSY that are beyond the scope of this review. 

For these reasons, this review is devoted to searches 
for SUSY at colliders. In this Section, we present brief 
summaries of the supersymmetric spectrum, parameters, 
and unifying frameworks to establish our conventions and 
notation. More extensive phenomenological reviews of 
SUSY may be found in Refs. flU S3 Hi S3 . 



A. Superpartners 

In this review, we focus our attention on the min- 
imal supersymmetric extension of the standard model 
(MSSM), the supersymmetric model with minimal field 
content. Bosonic superpartners are given names with the 
prefix "s-," and fermionic superpartners are denoted by 
the suffix "-ino." Squarks and sleptons are collectively 
known as "sfermions," and the entire group of superpart- 
ner particles are often called "sparticles." 

The particle content of the MSSM is in fact slightly 
more than a doubling of the SM particle content. This 
is because, in addition to introducing superpartners for 
all known particles, the MSSM requires two electroweak 
Higgs doublets. There are two reasons for this. First, 
in the SM, mass terms are generated for up- and down- 
type particles by Yukawa couplings to ip* and ip, respec- 
tively, where ip is the SM Higgs field. In SUSY, Yukawa 
couplings are generalized to terms in a superpotential, 
a function of superfields that contain both SM parti- 
cles and their superpartners, which generates the SM 



Yukawa couplings as well as all other terms related to 
these by SUSY. Complex-conjugated fields are not al- 
lowed in the superpotential, however. As a result, two 
separate Higgs fields, denoted H u and Hd, are required 
to generate masses through the superpotential terms 



W = X U H U QU + X d H d QD + X e H d LE 



(5) 



where Q, U, D, L, and E are the SU(2) quark doublet, 
up-type quark singlet, down- type quark singlet, lepton 
doublet, and lepton singlet superfields, respectively, and 
the A couplings are Yukawa couplings. Second, SUSY re- 
quires that the SM Higgs field have fermion partners, the 
Higgsinos. The introduction of these additional fcrmions 
charged under SM gauge groups ruins anomaly cancella- 
tion, making this theory mathematically untenable. The 
introduction of an additional Higgs doublet, with its ex- 
tra Higgsinos, restores anomaly cancellation. 

The MSSM Higgs boson sector therefore consists of 
eight degrees of freedom. As in the SM, three of these 
are eaten to make massive W and Z bosons, but five 
remain, which form four physical particles: 



MSSM Higgs Bosons (Spin 0) : h, H, A, F d 



(6) 



where h and H are the CP-even neutral Higgs bosons, 
with h lighter than H, A is the CP-odd neutral Higgs 
boson, and H ± is the charged Higgs boson. 

The remaining supersymmetric particle content of the 
MSSM is straightforward to determine and consists of 
the following states: 



Neutralinos (Spin 1/2) 
Charginos (Spin 1/2) 

Sleptons (Spin 0) 

Squarks (Spin 0) 

Gluinos (Spin 1/2) 



B,W°,H°,H° d 

U-L,R, CL,R, th,R. 
d>L,R, 5l,r, ^l.R 
9 ■ 



(7) 



Each SM chiral fermion has a (complex) scalar partner, 
denoted by the appropriate chirality subscript. The di- 
mensionless couplings of all of these particles are fixed 
by SUSY to be identical to those of their SM partners. 
Note, however, that, as described in the appropriate sec- 
tions below, the states in each line of Eq. ^ (except 
for the last one) may mix, and mass eigenstates are in 
general linear combinations of these gauge eigenstates. 

Finally, most analyses of SUSY include the supersym- 
metric partner of the graviton: 



Gravitino (Spin 3/2) : G 



(8) 



Although not technically required as a part of the MSSM, 
when SUSY is promoted to a local symmetry, it necessar- 
ily includes gravity, and the resulting supergravity theo- 
ries include both gravitons and gravitinos. The gravitino 
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is therefore present if SUSY plays a role in unifying the 
SM with gravity, as in string theory. 

If SUSY were exact, the gravitino's properties would be 
determined precisely by the graviton's, and it would be 
massless and have gravitational couplings suppressed by 
the reduced Planck mass M* ~ 2.4 x 10 18 GeV. However, 
just as Goldstone bosons appear when conventional sym- 
metries are spontaneously broken, a fermion, the Gold- 
stino G1/2, appears when SUSY is broken. The gravitino 
then becomes massive by eating the Goldstino. In terms 
of F, the mass dimension-2 order parameter of SUSY 
breaking, the gravitino mass becomes 



F 



(9) 



and, very roughly, its interactions in processes probing 
energy scale E may be characterized by a dimensionless 
coupling 



96 



E 2 



E 2 



(10) 



Light gravitinos couple more strongly. As we will see 
below, in well-motivated supersymmetric theories, these 
properties may take values in the range 

eV < m d < 10 TeV (11) 

10" 5 > cjg >10" 18 , (12) 

where we have assumed colliders probing E ~ m wea k. 

B. Supersymmetry Parameters 

As noted above, if SUSY exists in nature, it must be 
broken. Although many different Lagrangian terms could 
be added to break SUSY, only some of these are allowed 
if SUSY is to stabilize the gauge hierarchy. These terms, 
known as "soft" SUSY-breaking terms, include most, but 
not all, Lagrangian terms with mass dimension 3 and 
below (48 ). For the MSSM, they are 



Q 



IQI 2 



-m%\D\ 2 - 



-m\\L\ 2 - 



M X BB + M 2 W 3 W j + M 3 g k g k 



m%\E\ 2 



h.c. 



\ u AuH u QU + \ d A D H d QD + X e A E H d LE 



\H d \ 2 + {BH u H d + h.c). (13) 



These lines are sfermion masses, gaugino masses, trilinear 
scalar couplings ( M-terms" ) , and Higgs boson couplings. 
In addition to the parameters above, there are two other 
key parameters: the fi parameter, which enters in the 
Higgsino mass terms \iH % u H d , and 



tan (3 = 



\ u 



(14) 



which parameterizes how the SM Higgs vacuum expec- 
tation value is divided between the two neutral Higgs 
scalars. 



The interactions of Eq. ( p~3| ) conserve i?-parity (|4U1 ISUjl . 
With R = (-i)3( s - i )+ 2 5' j where B and L are the baryon 
and lepton numbers, respectively, and S the spin, all su- 
perpartners are odd and all SM particles are even under 
i?-parity. This implies that all interactions involve an 
even number of superpartners, and so the lightest super- 
partner is stable, and a potential dark matter candidate. 
i?-parity violation generically violates both baryon and 
lepton number, leading to too-rapid proton decay, which 
is why, for most of this review, we limit ourselves to the 
i?-parity conserving case. 

Even restricting ourselves to the i?-parity-preserving 
terms of Eq. (13 1, however, we see that SUSY intro- 
duces many new parameters. Note that the terms in- 
volving sfermions need not be flavor-diagonal, and so the 
sfermion masses and A-terms are in fact matrices of pa- 
rameters in the most general case. At the same time, 
fully general flavor mixing terms violate low energy con- 
straints on flavor-changing neutral currents. In addition, 
arbitrary complex parameters also violate bounds on CP- 
violation from, for example, ex and the electric dipolc 
moments of the electron and neutron. These consider- 
ations motivate unifying frameworks, to which we now 
turn. 



C. Unifying Frameworks 

In collider searches, it is desirable to consider theo- 
ries that are both viable and simple enough to be ex- 
plored fully. For this reason, it is common to work 
in simple model frameworks that reduce the number of 
independent SUSY parameters. In some cases, these 
model frameworks also motivate particular collider sig- 
natures that might otherwise appear highly unlikely or 
fine-tuned. 

In the most common unifying frameworks, SUSY is as- 
sumed to be broken in some other sector. SUSY break- 
ing is then mediated to the MSSM through a mechanism 
that defines the framework. This sets SUSY-breaking 
parameters at some high energy scale. Renormalization 
group evolution to the weak scale then determines the 
physical soft SUSY-breaking parameters and the phys- 
ical spectrum of the MSSM. A representative example 
of renormalization group evolution is shown in Fig. [T] 
In this evolution from the high scale to the weak scale, 
gauge couplings increase masses and Yukawa couplings 
decrease masses. This is central to understanding the 
sparticle spectrum of many models. In addition, it ex- 
plains why m 2 H becomes negative at the weak scale — it 
is the only particle to receive large negative contributions 
from Yukawa couplings without compensating large pos- 
itive contributions from the strong coupling. When H u 
becomes tachyonic, it breaks electroweak symmetry, and 
this feature, known as "radiative electroweak symmetry 
breaking," is a virtue of many supersymmetric frame- 
works. Note, however, that radiative electroweak sym- 
metry breaking makes essential use of the large top quark 
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sleptoTis^- rn 
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log 10 (Q/1 GeV) 

FIG. 1 Renormalization group evolution of scalar and gaug- 
ino mass parameters from the GUT scale A/gut — 2 x 
10 16 GeV to the weak scale in a representative mSUGRA 
model. From Ref. (HU). 



mass, and so shifts the burden of understanding why elec- 
troweak symmetry is broken to the question of why the 
top quark is heavy. 

In this section, we discuss several common unifying 
frameworks that have been used in collider searches, 
namely, models with gravity-, gauge-, and anomaly- 
mediated SUSY breaking. Each of these has its dis- 
tinctive characteristics. As a rough guide, in Fig. [2] we 
show representative spectra resulting from each of these 
frameworks. These spectra may be generated using pub- 
licly available computer programs, including is A jet (jSTJ), 
SOFTSUSY flS2), SPHENO and SUSPECT (|54"|) . 



1. Gravity Mediation (SUGRA) 

In gravity-mediated SUSY-breaking models (fS"6l IST1 
l58| [59] l60l l6T]l . sometimes referred to as supergravity 
(SUGRA) models, SUSY breaking in a hidden sector is 
mediated to the MSSM through terms suppressed by the 
reduced Planck mass M*. For example, sfermion masses 
are ~ F/M*. For these to be at the weak scale, yF 
must be around 10 11 GeV. Given Eqs. ^ and ( fit)} , the 
gravitino also has a weak scale mass and couples with 
gravitational strength in SUGRA models. 

Without a quantum theory of gravity, the structure 
of gravity-mediated SUSY parameters is unconstrained 
and generically violates low-energy constraints. To make 
these theories viable, ad hoc unifying assumptions must 
be made. By far the most common assumptions are those 
of minimal supergravity (mSUGRA), which is specified 
by 4 continuous and 1 discrete parameter choice: 



mSUGRA: m , Jn 1 y 2 , A Ql tan/3, sign(/z) , 



(15) 



where the first three parameters are the universal scalar 
mass (including the two Higgs scalars), unified gaugino 
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FIG. 2 Sparticle spectra for representative models with 
gravity-, gauge-, and anomaly-mediated SUSY breaking. The 
masses are normalized to M2, the Wino mass parameter at the 
weak scale. In the gravity-mediated case, two example spec- 
tra are presented: one for "no-scale" models with mo = 0, 
and another for mo = 5M2. In the anomaly-mediated case, 
the sleptons are tachyonic in the minimal case — additional 
effects are required to raise these to a viable range. From 
Ref. (55J|. 



mass, and universal A-parameter, which are all speci- 
fied at the grand unified theory (GUT) scale Mqtjt — 
2 x 10 16 GeV. The remaining SUSY parameters \n\ and 
the dimension-2 Higgs boson mass parameter B are de- 
termined by requiring that the Higgs potential at the 
weak scale give correct electroweak symmetry breaking. 
At tree-level, this requires 



1 2 

2 z 

sin 2/3 



tan 2 /3 



tan 2 13 - 1 



2B 



H, 



+ m 



(16) 
(17) 



Gaugino mass unification is motivated by the unifica- 
tion of gauge couplings at Mqut in the MSSM. It leads 
to the prediction that the Bino, Wino, and gluino masses 
arc in the ratio Mi : M2 : M3 ~ 1 : 2 : 7 at the weak 
scale, as evident in Fig. [2] Scalar mass universality is on 
much less solid ground. Even in GUTs, for example, the 
Higgs scalars are not necessarily in the same multiplet 
as the squarks and sleptons. This motivates a slightly 
less restrictive framework, the non-universal Higgs model 
(NUHM) in which tuq is the universal sfermion mass, 
but n%H u and mji d are treated as independent parame- 
ters. One may exchange these new degrees of freedom 
for the more phenomenological parameters [i and tua at 
the weak scale: 



NUHM: m , m 1 / 2 , A , tan/3, fi, iua , 



(18) 



The NUHM framework is employed in some MSSM Higgs 
boson studies discussed in Sec. III.DI 
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2. GMSB 

In gauge-mediated SUSY-breaking (GMSB) mod- 
els ([62] |63l [64] E2 H2 E3, in addition to the gravity- 
mediated contributions to soft parameters discussed 
above, each sparticle receives contributions to its mass 
determined by its gauge quantum numbers. These 
new contributions to sfermion masses are ~ F/M mess , 
where M mcss is the mass scale of the messenger parti- 
cles that transmit the SUSY breaking. The GMSB con- 
tributions are flavor-blind, and do not violate low en- 
ergy bounds. For these to be dominant, one requires 
M mess < 10 14 GeV, and so we find that m & ~ F/M* « 
F/M m e S s ~ m-wcak in GMSB scenarios, and the lightest 
supersymmetric particle (LSP) is always the gravitino. 

In GMSB models, the collider signatures are deter- 
mined by the next-to-lightest supersymmetric particle 
(NLSP) and its lifetime, or equivalently, the gravitino 's 
mass. If the NLSP is the lightest neutralino, the col- 
lider signature is either missing energy or prompt pho- 
tons, Z or Higgs bosons from x° — > (7i Z, h)G ((551 159) 1: 
if the NLSP is a slepton, the signature is typically ei- 
ther long-lived heavy charged particles or multi-lepton 
events (i70| l7T| l72| . 



3. AMSB 

A third class of SUSY models are those with anomaly- 
mediated SUSY-breaking (AMSB) (|73] [BJ). These are 
extra dimensional scenarios in which SUSY is broken 
on another 3-dimensional subspace, and transmitted to 
our world through the conformal anomaly. As with all 
anomalies, this effect is one-loop suppressed. The funda- 
mental scale of SUSY breaking as characterized by the 
gravitino mass is therefore mQ ~ 10 — 100 TeV, with 
MSSM sparticle masses one-loop suppressed and at the 
weak scale. 



The AMSB contributions to sparticle masses are com- 
pletely determined by the sparticle 's gauge and Yukawa 
couplings. This leads to a highly predictive spectrum. 
Unfortunately, one of these predictions is rr?z,rr?~ < 0, 
but various mechanisms have been proposed to solve this 
tachyonic slepton problem; see, e.g., Refs. (|75l l76l l77)l . 

The gaugino masses are determined by the correspond- 
ing gauge group beta functions. In particular, AMSB 
predicts M% : M 2 : M 3 ~ 2.8 : 1 : 8; because the 
SU(2) coupling is nearly scale-invariant in the MSSM, 
the Wino mass is the smallest. AMSB scenarios there- 
fore motivate supersymmetric models with W LSP and 
W ± NLSP. This triplet may be extremely degenerate, 
with the chargino traveling macroscopic distances before 
decaying to soft and invisible decay products, which pro- 
vides a distinctive and challenging signature for collider 
searches (fT8"jl . 



D. Supersymmetric Higgs Bosons 

The MSSM Higgs potential is 

Vh = (mi u + | A i| a )|fl2| a + (mi ri + |A t | a )|flS| a 

- sCiM + hxO + ^W^T 



1 



+ Ug* + ^)(\HW-\H«\*Y 



(19) 



where the parameters [i, 171 'h u j and 771 lf d are as defined in 
Sec. |II.B| and SUSY implies that all the quartic couplings 
are determined by the SU(2) and U(l) hypercharge gauge 
couplings, denoted g and g', respectively. Vh automati- 
cally conserves CP, since any phase in the B parameter 
can be eliminated by a redefinition of the Higgs fields. 

Assuming these parameters are such that the potential 
admits a stable, symmetry-breaking minimum, the three 



parameter combinations m H 



|/i| 2 , and B 



may be exchanged for the two vacuum expectation val- 
ues v u = \/2{H®), v d = v2(-ff°), and one physical Higgs 
boson mass, conveniently taken to be rriA- The W bo- 
son mass fixes v 2 + v 2 ,, leaving one additional degree of 
freedom, usually taken to be tan/3 = v u /vd- Thus, at 
tree-level, the entire MSSM Higgs boson sector is deter- 
mined by two parameters, tua and tan (3. 

In terms of these parameters, the physical Higgs boson 
masses are 



m\ + 171% i ^J(m\+m%)' 2 — Am A m 2 z c^p 



2 2,2 

m H± = m A + m w , 



(20) 
(21) 



where C2/3 = cos 2(3. The CP-even mass eigenstates are 
related to the gauge eigenstates through 

H\_f cos a sina\ / y/2ReH% - v d \ , . 
h J \ -sina cosa J \ \/2ReH° - v u J ' 1 ' 

where the rotation angle a satisfies 



cos 2a = — cos 2(3- 



1 A ,,l Z 



(23) 



with — tt/2 < a < 0. 

Equation ( p0| ) implies that m,h < mz \ cos2/3|, a rather 
disastrous relation, given that experimental bounds ex- 
clude rrih < mz- The results presented so far, however, 
are valid only at tree-level. Large radiative corrections 
from top squark/quark loops (see, e.g., Ref. (79)), 



A? 



1 



3g 2 mj 



1 (3 8tt 2 ? 



log 



(24) 



can lift TUh to values above the experimental bounds. 
Note, however, that for tan/3 = 1, rrih — at tree-level, 
and so large values of rrih are not possible for tan [3 ~ 1 . 
From considerations of the Higgs mass alone, tan/3 < 1 
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is possible. However, such values imply very large top 
Yukawa couplings, which become infinite well below the 
GUT or Planck scales. In addition, in simple frameworks, 
tan (3 < 1 is incompatible with radiative electroweak sym- 
metry breaking; for a review of bounds on tan/3, see 
Ref. ([Mi 



E. Neutratinos and Charginos 

The neutralinos and charginos of the MSSM are the 
mass cigcnstatcs that result from the mixing of the elec- 
troweak gauginos B and with the Higgsinos. 

The neutral mass terms are 



1 



{^) T M N ^ + h.c. 



(25) 



where (i/>°) T = {-iB, -iW 3 , H% H°) and 
/ Mi 



M 



N 



Mi 


\g'vd 



M 2 
\gvd 
-\gv u 



-\g'vd 
\gvd 




2 



g'v u \ 
\gv u 

Q 



The neutralino mass eigenstates are Xi — ^ijWji where 
N diagonalizes Mjv- In order of increasing mass, the four 
neutralinos are labeled Xi, xSj; an d Xa- 
The charged mass terms are 



(V~) T M c V + + h.c. , 
where {ij) ± ) T = (-iW ± ,H ± ) and 



(27) 



These masses are given in third-generation notation; in 
the presence of flavor mixing, these generalize to full six- 
by-six and three-by-three matrices. 

The left-right mixing is proportional to lepton mass, 
and is therefore expected to be insignificant for sclcctrons 
and smuons, but may be important for staus, especially if 
tan/3 is large. Through level repulsion, this mixing low- 
ers the lighter stau's mass. As noted in Sec. |II.C] Yukawa 
couplings also lower scalar masses through renormaliza- 
tion group evolution. Both of these effects imply that in 
many scenarios, the lighter stau is the lightest slepton, 
and often the lightest sfcrmion. 



G. Squarks 



The mass matrix for up-type squarks is 
/m^+mf+m|(|-fs^)c 2/ 3 m t (A t -(J,cot j3) 



m t (A t - /j, cot /3) 



(26) in the basis (iz,,tii), and for down-type squarks is 



(31) 



/m|+mf-mf (|-|sfj,)c 2/ a m b (A b -fita,n(3) \ 
I m b (A b -fi tan (3) m 2 £) +ml-m%\s^ v C2pJ 

( 32 ) 

in the basis (6l,^_r)- Large mixing is expected in the 
stop sector, and possibly also in the sbottom sector if 
tan/3 is large. Because of these mixings and the impact 
of large Yukawa couplings in renormalization group evo- 
lution, the 3rd generation squarks are the lightest squarks 
in many models (fBTjl . 



M, 



M, 



c 



V2 



gv u 



V2 



gvd 



(28) 



Vijip^ and 



The chargino mass eigenstates are Xi — y 13 <Vj 
Xl = \Jijip~, where the unitary matrices U and V are 
chosen to diagonalize Mc, and xt is lighter than xf- 

F. Sleptons 

Sleptons are promising targets for colliders, as they are 
among the lightest sparticles in many models. As noted 
in Sec. |II.A[ sleptons include both left- and right-handed 
charged sleptons and sneutrinos. The mass matrix for 
the charged sleptons is 

t|+ra2_ TO |(i _ s 2^) C2/3 m r (yl r -/itan/3) 
m T (A T -/xtan/3) m|+rn2_ TO | s 2^ C2/3 

(29) 

in the basis (t£,77j), where sw = sin#w The sneutrino 
has mass 



m| = m| + -m|cos2/3 . 



(30) 



III. SEARCHES FOR MSSM NEUTRAL HIGGS BOSONS 

As already explained in Sec. [IT] two Higgs doublets are 
needed in the MSSM to give mass to both up- and down- 
type quarks. Under the assumption that the Higgs sec- 
tor is CP conserving, the physical states are two neutral 
CP-even Higgs bosons (h and H, ordered by increasing 
mass), a neutral CP-odd Higgs boson (A), and a dou- 
blet of charged Higgs bosons (H ). Further details on 
the Higgs sector of the MSSM have been given in Sec. [H] 
Here, we focus on searches for the neutral Higgs bosons 
of the MSSM, while searches for charged Higgs bosons 
will be discussed in Sec. HVl 



A. MSSM benchmark scenarios 

It has been seen that two parameters are sufficient 
to fully determine the MSSM Higgs sector at tree level. 
These are commonly taken to be the A boson mass tua 
and tan /3, the ratio of the vacuum expectation values of 
the Higgs fields giving mass to the up- and down-type 
quarks. This picture is modified significantly, however, 
by large radiative corrections, arising essentially from an 
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incomplete cancellation of the top and stop loops. In 
particular, the important prediction < mz\cos2(3\ is 
invalidated. Among the many parameters of the MSSM, 
a few have been identified as being most relevant for the 
determination of Higgs boson properties. In addition to 
rriA and tan f3, an effective SUSY breaking scalar mass, 
MgjjsY, which sets the scale of all squark masses, and a 
term controlling the amount of mixing in the stop sector, 
X t , play the leading role. (In Eq. (24 1, the stop mass is 



directly related to Msusy, and stop mixing is neglected.) 
The model is further specified by a weak gaugino mass, 
M2, the gluino mass, nig, and the SUSY Higgs mass term 
/i. The relation X t = A — /icot/3 then allows the trilin- 
ear Higgs-squark coupling A (assumed to be universal) 
to be calculated. For large values of tan/3, mixing in the 
sbottom sector becomes relevant too; it is controlled by 
Xi, = A — (i tan [3. Finally, the top quark mass m t needs 
to be specified. 

A few benchmark scenarios (82; 83) were agreed upon 
to interpret the searches for MSSM Higgs bosons. The 
most widely considered are the so-called "m/j-max" and 
"no-mixing" ones, where Mgugy = 1 TeV, M 2 — 
200 GeV, n = -200 GeV and m § = 800 GeV. In m h - 
max, X t is set equal to 2M$usy (in the on-shell renor- 
malization scheme), while it is set to in the no-mixing 
scenario. The largest value of is obtained for large 
rriA and tan (3, and is maximized (minimized) in the m^- 
max (no-mixing) scenario. In the m/j-max scenario, the 
maximum value of is ~ 135 GeV. 



B. Searches at LEP 

At LEP, the neutral Higgs bosons of the MSSM have 
been searched for in two production processes, the Hig- 



gsstrahlung process e 



hZ (JHU |85l [86]) and the 
associated production e + e~ — » hA (j8~7J). Both processes 
are mediated by s-channel Z boson exchange. With the 
notations of Sec. ITT1 the cross sections are 



<i hz = sin 2 {(3- a)alf 
<r hA = cos 2 (/3 - a)\c7hz 



where /3 and a are defined in Eqs. (14) and (22 1, 



SM 
>hZ 



If,, 2 + n 2 U i/2 *hz + 12m|/a 
967TS 1 e + e) hz (1 - ml/3) 2 



(33) 
(34) 



(35) 



is the SM Higgs boson production cross section, s is the 
square of the center-of-mass energy, and 

Kj = [1 - (m + mj) 2 /s][l - (mi - rrij) 2 / s] , (36) 



(37) 



It is apparent from the above formulae that the two 
processes are complementary. In practice, the Hig- 
gsstrahlung process dominates for values of tan (3 close 
to unity, while associated production dominates for large 
values of tan (3, if kinematically allowed. In large regions 



of the MSSM parameter space, the h decay branching 
fractions are similar to those of the SM Higgs boson. For 
a mass of 115 GeV, these are 74% into bb, 7% into both 
t+t- and gg, 8% into WW*, and 4% into cc (JHHJ) ■ The A 
boson couples only to fermions, so that its decay branch- 
ing fraction into bb is always close to 90%, with most of 
the rest going into t + t~ . These same branching fractions 
also hold for the h boson for large values of tan (3 (|55|) . 

Searches for Higgs bosons were performed at LEP first 
in Z boson decays during the LEP1 era, and subse- 
quently at increasing center-of-mass energies at LEP2, up 
to 209 GeV in 2000. In the following, only the searches 
performed at the highest energies are described. 

The four LEP experiments carried out searches for the 
SM Higgs boson produced via Higgsstrahlung, e + e~ — > 
HZ, and the results were combined to maximize the sen- 
sitivity. 1 Four final state topologies were analyzed to 
cope with the various decay modes of the Higgs and 
Z bosons: a four-jet topology with two 6-tagged jets, 
for (H — > bb)(Z — > qq); a two fr-tagged jets and two- 
lepton topology, for (H — > bb)(Z — > £ + t~), with £ = e or 
fi; a two 6-tagged jets and missing energy topology, for 
(H — > bb)(Z — > vv)\ and a two-jet and two-r topology 
for (H -> bb)(Z -> t+t-) and (H -> t+ T -){Z -> qq). 
A few candidate events were observed at the edge of the 
sensitivity domain, but the overall significance was only 
at the level of 1.7 a. A lower mass limit was therefore 
derived, excluding a SM Higgs boson with mass smaller 
than 114.4 GeV (|90j). 

The Higgs boson mass lower limit depends on the 
strength of the HZZ coupling, and the LEP collabo- 
rations also provided, as a function of the mass of a 
SM-like Higgs boson, an upper limit on £ 2 , where £ is 
a multiplicative factor by which the SM HZZ coupling 
is reduced (90). By SM-like, it is meant that the decay 
branching fractions are similar to those expected from a 
SM Higgs boson. This result is shown in Fig. [3] Con- 
straints on the MSSM parameter space can be deduced 
from this, since in that case £ = sin(/3 — a). 

For Higgs boson masses accessible at LEP, the struc- 
ture of the MSSM Higgs sector is such that the h and 
A masses are similar whenever associated production is 
relevant, i.e., for large values of tan/3. Searches for hA 
associated production were performed in the four 6-jet 
final state for (h — > bb) (A — > bb) and in the two 6-jet and 
two-r topology for (h/A — > bb)(A/h — > t + t~). The con- 
straint that the h and A boson candidate masses should 
be similar was imposed. The backgrounds from multijet 
and WW production were largely reduced by the fr-jet 
identification requirements, leaving ZZ as an irreducible 
background. 

No significant excess over the SM background expec- 



1 Production by vector boson fusion, e + e~ — > He + e~ or 
Hvv H89j l. was also considered, but its contribution was found 
to be negligible in practice. 
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120 

m H (GeV/c 2 ) 

FIG. 3 The upper bound on the factor £ 2 by which the square 
of the SM HZZ coupling is multiplied, as provided by the 
LEP experiments l|90p. The full curve is the observed limit, 
the dashed curve the median expected limit in the absence of 
signal, and the green and yellow bands are the 68% and 95% 
probability regions around the expected limit. 



tation was observed, and production cross section upper 
limits were derived as a function of to,/, ~ m^. For each 
benchmark scenario, a scan was performed as a function 
of tyia and tan [3, and in each point of the scan the cross 
section upper limit was compared to the corresponding 
prediction, taking into account the slight modifications 
expected for the values of the h and A branching fractions 
into bb and t + t~, as well as the non-negligible difference 
between and which develops at lower values of 
tan (3. If the cross section upper limit was found to be 
smaller than the prediction, the (m a, taxi [3) set was de- 
clared excluded. The result of the combination of the 
searches in the hZ and hA channels by the four LEP 
experiments ([51)1 is shown in Fig. [1] projected onto the 
(to/,, tan/3) plane in the m/,-max and no- mixing scenar- 
ios. In the derivation of those results, contributions of 
the e + e~ — > HZ and HA processes were also taken into 
account whenever relevant, where H is the heavier CP- 
even Higgs boson. 

In the most conservative scenario, i.e., m^-max, it can 
be seen in Fig. [4] that the lower limit on the mass of the 
SM Higgs boson holds also for to/, as long as tan/3 is 
smaller than about 5, and that values of tan/3 between 
~ 0.7 and 2 are excluded for the current average value of 
the top quark mass, 173.1 ±1.3 GeV (92). A lower mass 
limit of 93 GeV is obtained for to/, ~ vha for large values 
of tan /?. 

The benchmark scenarios were chosen such that the 
Higgs bosons do not decay into SUSY particles. An in- 
teresting possibility is that the h — > XiXi decay mode is 
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FIG. 4 Domains excluded at 95% CL (light green) and 99.7% 
CL (dark green) by the four LEP experiments (|91f) in the 
(m/,,tan,3) plane in the mj-max (top) and no-mixing (bot- 
tom) benchmark scenarios, with mt = 174.3 GeV. The yellow 
regions are not accessible theoretically. The dashed lines rep- 
resent the boundaries of the domains expected to be excluded 
at 95% CL in the absence of signal. The upper boundaries of 
the physical regions are indicated for four top quark masses: 
169.3, 174.3, 179.3, and 183 GeV, from left to right. 



kinematically allowed, where x? is the LSP. If i?-parity 
is conserved, the LSP is stable and, since it is weakly 
interacting, the Higgs boson decay final state is invisi- 
ble. Searches for such an "invisible" Higgs boson were 
performed by the LEP experiments, and the combina- 
tion (|93|l yields a mass lower limit identical to that set 
on the SM Higgs boson if the production cross section is 
the SM one, as is the case for low values of tan j3. 
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To cope with fine-tuned choices of MSSM parameters, 
the LEP collaborations considered yet other possibilities, 
e.g., that the h — ► AA decay mode is kinematically al- 
lowed, or that the h — > bb decay is suppressed. For 
example, dedicated searches for hA — > AAA — > bbbbbb 
and for ft,Z, with /i -> in a flavor-independent way, 
have been performed (jMJ). In the end, the sensitivity of 
the standard searches is only slightly reduced, except for 
rather extreme parameter choices leading, for instance, to 
nih — 100 GeV, while at the same time tua < 2m&. This 
last possibility is however less unnatural in extensions 
of the MSSM, such as the NMSSM where an additional 
Higgs singlet field is introduced (j9"S"|) . 

Finally, the possibility that CP is violated in the Higgs 
sector has also been considered. While CP is conserved 
at tree level, radiative corrections may introduce such 
a CP violation if the relative phase of /i and A is not 
vanishing. In such a case, the three mass eigenstates all 
share properties of h, H and A, so that the signatures of 
Higgs boson production are less distinct. The constraints 
are accordingly weaker. A dedicated "CPX" scenario (j96l 
I97p was set up to perform quantitative studies. As an 
example, a region around rrih = 45 GeV and tan /3 = 5 
is not excluded for M SUSY = 500 GeV, M 2 = 200 GeV, 
H = 2 TeV, and m s = 1 TeV, when \A\ = 1 TeV and 
arg(A) = 90°. Further details can be found in Ref. ([91). 



C. Searches at the Tevatron 

At the Tevatron, i.e., in pp collisions at 1.96 TeV, the 
dominant production mechanism for the SM Higgs bo- 
son is via gluon fusion, gg — > H (l98| l99|l . In the mass 
range that is of interest for a SM-like Higgs boson of 
the MSSM, namely rrih < 135 GeV, the dominant de- 
cay mode is H — > bb. Such a two-jet final state is to- 
tally overwhelmed by standard jet production via the 
strong interaction, even after 5-jet identification. This is 
why the SM Higgs boson searches at the Tevatron have 
been performed in the associated production processes 
qq — ► (W/Z)H (llOOp . which proceed via s-channel W 
or Z exchanges in a similar way to the Higgsstrahlung 
in e + e~ collisions. In spite of cross sections an order 
of magnitude smaller than that of gluon fusion, these 
processes offer better discrimination against the multijet 
background, by making use of the leptonic decays of the 
W and Z (W -► tv, Z -> t+f and Z -> w). These 
searches for the SM Higgs boson apply equally well for 
the h boson of the MSSM in the low tan/3 regime. Their 
sensitivity is, however, still not sufficient to provide any 
significant constraint. 

The situation is much more favorable for large values of 
tan f3. In this regime, the A boson is almost mass degen- 
erate with cither the h or H boson, depending on whether 
iriA is less than or greater than m™ ax , where m™ ax is the 
maximum value that rrih can take, e.g., 135 GeV in the 
m/j-max scenario. In the following, the two nearly de- 
generate Higgs bosons are collectively denoted 4>. Their 



couplings to b quarks and r leptons are enhanced by a fac- 
tor tan/3 with respect to the SM couplings. As a result, 
the contribution of the b quark loop to their production 
via gluon fusion is enhanced by a factor 2 tan 2 /3. Al- 
though this is not sufficient to render feasible a detection 
in the </> — > bb decay mode, this is not the case for the 
4> — > t + t~ decay mode, which has a branching fraction 
of ~ 10%. 

Both CDF and D0 required one of the two r leptons 
to decay leptonically (r — > [ej [i)vv) to ensure proper 
triggering. Three final state topologies were considered: 
e^had, M r had, and efj,, all with missing transverse energy 
IfjT from the r decay neutrinos. Here Thad denotes a r 
lepton decaying into hadrons and a neutrino. The dom- 
inant, irreducible background comes from Z production 
with Z — ► t + t~, but there also remains a substantial 
component from (W — > £v)+jet, where the jet is misiden- 
tified as a r lepton. This background was reduced, for 
instance, by requiring a low transverse mass of the lepton 
and the $t ■ The final discriminating variable was cho- 
sen to be the visible mass m V i S = \J (P Tl + Pt 2 +Pt) 2 , 
constructed from the r visible products and from the 
IpT- The distribution of m v i s obtained by CDF (jlOlf) in 
a 1.8 fb _1 data sample is shown in Fig. [5} From this dis- 
tribution, as well as from a similar one in the e^i channel, 
a cross section upper limit on <fi production was derived, 
which in turn was translated into exclusion domains in 
the (mA, tan (3) plane within benchmark scenarios. The 
result obtained in the m/i-max and no-mixing scenarios 
is shown in Fig. [6] Similar results have been obtained by 
D0 {TD5J. The calculations of Ref. (TTM]l were used to 
derive these results as well as those reported in the rest 
of this section. 

Because of the enhanced coupling of <f> to b quarks at 
high tan /3, the production of Higgs bosons radiated off a 
b quark may be detectable in the <j> — > bb decay mode in 
spite of the large background from multijet events pro- 
duced via the strong interaction ("QCD background"). 
This process can be described in the so-called four-flavor 
or five-flavor schemes, and it has been shown that the two 
approaches yield very similar results (104). In the four- 
flavor scheme, the main contribution comes from gluon 
fusion, gg — > bb<p, while the main one in the five-flavor 
scheme comes from gb — > b(f>. Because one of the fi- 
nal state b quarks (a spectator b quark in the five-flavor 
scheme) tends to be emitted with a low transverse mo- 
mentum, the searches required only three b jets to be 
identified. The signal was searched for by inspecting the 
mass distribution of the two jets with highest transverse 
momenta in the sample of events with three 6-tagged jets. 
Further discrimination against the QCD background was 
provided by the mass of the charged particles in the 
tagged jets (at CDF (105)) or by the inclusion of ad- 
ditional kinematic variables in a likelihood discriminant 
(at D0 (fi"06"10 . The QCD background was modeled us- 
ing a combination of information from control samples in 
the data, where one of the jets is not 6-tagged, and from 
Monte Carlo simulations of the various processes con- 
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FIG. 5 Visible mass distribution in the (e//i)rhad channels 
from the CDF search for 4> — > TT (101). The signal contri- 
bution indicated corresponds to the cross section upper limit 
set with this data. 
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FIG. 6 Domains in the (m.A, tan /3) plane excluded by the 
CDF search for 4> — > tt (|101|) . The domains excluded at LEP 
are also indicated. 
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FIG. 7 Fit to the mass of the two jets with highest trans- 
verse momenta in the CDF sample of events with three b- 
tagged jets (|105[) . The contributions of the various multijct 
backgrounds and of a signal with a mass of 150 GeV are in- 
dicated. 



tributing to the background (bbb, bbc, bbq, ccc, ccq, etc., 
where q represents a light quark, u, d, s, or a gluon). 
The mass distribution obtained by CDF in a 1.9 fb -1 
data sample is shown in Fig. [7] with the individual back- 
ground contributions displayed. No signal was observed, 
and production cross section upper limits were derived, 
from which exclusion domains in the (m^tan/?) plane 
were determined in various benchmark scenarios. The 
D0 result obtained with 2.6 fb -1 of data in the m^- 
max scenario is shown in Fig(8j In the derivation of the 
cross section upper limits and exclusion domains, special 
attention was given to a proper handling of the Higgs 
boson width, which is enhanced by a factor tan 2 /? at 
tree level and therefore becomes large with respect to 
the mass resolution. (This is not the case for <p ^ t + t~ 
because of the degradation of the mass resolution due to 
the missing neutrinos.). It should also be noted that the 
exclusion domain is quite sensitive to the model parame- 
ters. It is smaller in the no-mixing scenario, and also if fi 
is positive. These effects due to potentially large SUSY 
loop corrections to the production cross sections and de- 
cay widths tend to cancel in the search for tfi — > t + t~ 
described above (jWl ITUH]) . 

Finally, Higgs bosons produced in association with b 
quarks can also be searched for in the (f> — ► t + t~ decay 
mode. Although the branching fraction is an order of 
magnitude smaller than the one of (f> — > bb, the signal is 
much easier to disentangle from the background. A D0 
analysis (|109p was performed where one of the r leptons 
decays into a muon and neutrinos, while the other decays 
into hadrons and a neutrino. Furthermore, a 6-tagged 
jet was required, at which point the main background 
comes from top quark pair production, tt — > [ivbrvb. A 
neural network was used to discriminate signal and tt 
background, taking advantage of the large differences in 
their kinematic properties. The result, based on 1.2 fb -1 
of data, is shown in Fig. [9] Given the limited amount 
of integrated luminosity used up to now, this channel 
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FIG. 8 Domain in the (mA, tan /?) plane excluded by the 
D0 search for <j> — *• bb in events with three b-tagged jets in 
the m^-max scenario (|106p . 
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FIG. 9 Domain in the (mA, tan /3) plane excluded by the D0 
search for 4> —> t + t~ produced in association with a 6 quark 
in the mk-max scenario (j 109ft . The light green shaded region 
is the extension of the LEP exclusion region to tan/? > 50. 



appears to be quite promising. 



IV. SEARCHES FOR CHARGED HIGGS BOSONS 

Many extensions of the SM involve more than one com- 
plex doublet of Higgs fields. Two-Higgs doublet mod- 
els (2HDMs) fall into three main categories. In Type I 
models, all quarks and leptons couple to the same Higgs 
doublet. In Type II models, down-type fermions couple 



to the first Higgs doublet, and up- type fermions couple 
to the second Higgs doublet. Flavor-changing neutral 
currents are naturally avoided in Type I and Type II 
2HDMs. In Type III models, fermions couple to both 
doublets, and flavor-changing neutral currents must be 
avoided using other strategies. In addition to the three 
neutral Higgs bosons discussed in the previous section, 
2HDMs involve a pair of charged Higgs bosons, H ± . 
Most of the experimental results on charged Higgs bosons 
have been obtained within the context of Type II 2HDMs, 
of which the MSSM is a specific instance. Further de- 
tails on extended Higgs boson sectors may be found in 
Ref. ffTTTn ITTTTl . 

In Type II 2HDMs, the charged Higgs boson decay 
width into a fermion pair f u fd is 

x (m 2 fd tan 2 (i + m 2 fu cot 2 (3) , (38) 

where N c is the number of colors, and we have ap- 
proximated mf d <C m H ± in the phase space factor. 
Charged Higgs bosons therefore decay into the heaviest 
kinematically-allowed fermions: t~v t at large tan/3 and 
cs at low tan (3 for charged Higgs boson masses to which 
current accelerators are sensitive. 



A. Searches at LEP 

At LEP, charged Higgs bosons are produced in pairs 
through e + ' eT — > H + H~ (|112p . The production cross 
section depends only on SM parameters and on the mass 
of the charged Higgs boson. The process e + er — > H + W~ 
has a significantly lower cross section. 

The charged Higgs boson can decay into cs or tv t . 
In searches for Type I 2HDM Higgs bosons, the de- 
cay — > AW ± * (I113p was also considered, as in 
Refs. (|114l I115"]l . The interpretation of the search results 
generally assumed that Br(H ± — > tv t ) + Br(iJ ± — > qq') 
= 1, where the dominant qq' flavors are cs, due to the 
Cabbibo suppression of cb. This assumption leads to 
the consideration of three topologies for pair-produced 
charged Higgs bosons: four jets from H + H~ — > cscs, two 
jets, a t lepton and missing energy from H + H~ — > cstv t 
and two charge conjugate, acoplanar 2 r leptons from 

H+H~ -> T + V T T~V T . 

Direct searches for pair production of charged Higgs 
bosons have been published by all four LEP experi- 
ments (|114l 11151 11161 1117)1 . Each topological analysis 
began with a general selection for the expected num- 
ber of jets and r leptons, followed by more sophisti- 



2 The acoplanarity angle is the angle between the projections of 
the t momenta on a plane transverse to the beam axis. If this 
angle is less than 180°, the r leptons are said to be acoplanar. 
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FIG. 11 For a charged Higgs boson mass ol 100 GeV, the 
branching ratios for the top quark decay into H + b, under the 
assumption that Br(t -* W + b)+Br(t -> H + b)=l, and for 
the various H + decay channels, as a function of tan/3. From 
Ref. ipofl . 



FIG. 10 Limit on the charged Higgs boson mass as a function 
of Br(H + — ^ t + v t ), from the combined data of the four LEP 
experiments at center-of-mass energies from 189 to 209 GeV. 
The expected exclusion limit is shown as a thin solid line and 
the observed limit as a thick solid line; the shaded region is 
excluded dTTHJ. 



cated techniques. The main difficulty in these analy- 
ses was separating the signal from the nearly identi- 
cal signature of W + W~ production; selection criteria 
usually included a mass-dependent optimization. Tech- 
niques such as linear discriminants, likelihood estimators, 
and jet-flavor tagging were used in these analyses. The 
H + H~ — > r + v t t~ v T channel had additional complexity 
due to the missing neutrinos, which removed the possi- 
bility of reconstructing the -ff candidate masses and of 
improved discrimination from the equal-mass constraint. 
However, final states with t leptons can benefit from ex- 
tracting information about their polarization; the r + lep- 
ton from a H + boson (a scalar) is produced in a helicity 
state opposite to that of a t + lepton from W + decay. 

The LEP experiments have combined the results of 
their searches for charged Higgs bosons into one re- 
sult (|118p based on common assumptions. The to- 
tal dataset has an integrated luminosity of 2.5 fb _1 , 
collected at center-of-mass energies between 189 and 
209 GeV. The possible decays were restricted to H + — ► cs 
and t + v t in a general 2HDM framework. The com- 
bined mass limit is shown in Fig. [10] as a function of 
Br(_ff + — > t + v t ). A lower bound of 78.6 GeV holds for 
any value of the branching ratio. 



B. Searches at the Tevatron 

At the Tevatron, pair production of charged Higgs 
bosons is expected to occur at a very low rate. How- 
ever, in contrast to searches at LEP, advantage can be 
taken of the large mass of the top quark, which opens new 
ways to search for evidence of charged Higgs bosons. Two 
approaches have been considered, depending on whether 
the charged Higgs boson is lighter or heavier than the 
top quark. In the first case, the top quark can decay into 
a H + boson and a b quark (|112[1 . For heavier charged 
Higgs bosons, resonant production of a single H + boson 
followed by the decay H + — > tb is the most promising 
process f|119j) . 

In the SM, the top quark decays almost exclusively into 
a W boson and a b quark, and the possible signatures of 
tt pair production are associated with the various combi- 
nations of VF-boson decay channels. If the charged Higgs 
boson is lighter than the top quark, the decay t — > H + b 
will compete with the standard t — > W + b mode. The 
decay of the charged Higgs boson, with branching ra- 
tios different from those of the W boson, will modify 
the fractions of events observed in the various topolo- 
gies, compared to the SM expectations. 3 The qualita- 
tive aspects and magnitude of these modifications depend 
on the model parameters. The dependence on tan/3 of 
the top quark decay ratio to H + b and of the various 
charged-Higgs boson decay channels is shown in Fig. [TT] 



3 The branching ratios for topologies arising from SM tt pair pro- 
duction are roughly 50% in six jets; 14% in each of e, fi, and 
t + 4 jets +IpT\ 1% i n each of ee, fi/i, and tt + 2 jets +IpT\ & n d 
2% in each of efi, er, and /it + 2 jets +IpT- In each of these 
channels, there are two b jets. 
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FIG. 12 Limit on the mass of the charged Higgs boson as 
a function of tan j3 from the D0 search in top quark de- 
cays ipoji . 



for m H + = 100 GeV and a typical set of MSSM param- 
eters, with QCD, SUSY-QCD and electroweak radiative 
corrections to the top and bottom quark Yukawa cou- 
plings calculated with the CPsuperH code (j!2ip . The 
dominant H+ boson decay channels arc cs at low values 
of tan (3 and t+v t at high values; with this set of param- 
eters, H+ boson decays to W+ A/h are also allowed, al- 
though always at a small rate. The H+ — > t*b — > W+bb 
decay mode becomes relevant for charged Higgs boson 
masses closer to the top quark mass. It can be seen 
that charged Higgs bosons will be most prominent at 
high and low values of tan (5. Two simplified models 
address each of these regions: the tauonic model, with 
Br(H+ — ► t+v t )=1, and the leptophobic model, with 
Br(H+ — > cs)=l. The tauonic model is a very good ap- 
proximation to the MSSM with tan/3 > 15, while purely 
leptophobic charged Higgs bosons can be found in some 
multi-Higgs-doublet models (|122p . 

Analyses based on measurements of ti final states in- 
clude an earlier CDF search in 200 pb _1 of data (|123p 
and a recent D0 analysis of 1 fb _1 of data (120). The 
yields observed in the various topologies were compared 
to what would be expected in models with charged Higgs 
bosons, taking into account the t — > H + b and H+ de- 
cay branching ratios predicted as a function of the Higgs 
boson mass and tan (3. In particular, no excess of final 
states involving r leptons was observed, nor was any dis- 
appearance of final states with one or two leptons, jets 
and^?x, as would be expected at large and small tan/3, 
respectively. Figure [12] displays the exclusion domain in 
the plane of the charged Higgs boson mass and tan (3 from 
the D0 analysis Q120p . for leptophobic and tauonic mod- 
els. The CDF analysis excludes Br(i -> H+b) > 0.4 for 
a tauonic H+ boson (|123|l . 

In a recent analysis based on a data sample of 2.2 fb -1 , 
the CDF collaboration used a different approach to 
search for a leptophobic charged Higgs boson in top 
quark decays (124). The search was performed in the 
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FIG. 13 For a leptophobic charged Higgs boson, upper limit 
on the branching ratio Br(t — * H+b) as a function of the Higgs 
boson mass from a CDF search in top quark decays (| 124[l . 



lepton+jets+^r final states with two 6-tagged jets, where 
the lepton (electron or muon), the neutrino (responsible 
for the missing Ex), and a & jet were the signature of a 
t — ► Wb — > ivb decay, while the other top quark of the 
ti pair was assumed to decay to either Wb — ► qq'b or 
Hb — > csb. The tt events were fully reconstructed, taking 
the masses of the W boson and of the top quark into 
account as constraints to assign correctly each of the b 
jets to its parent t or t. Templates of the mass of the 
dijet system reconstructed from the non-& jets were used 
to extract limits on the branching ratio of t — > H+b, as 



shown in Fig. 13 



If the charged Higgs boson is heavier than the top 
quark, it will decay dominantly into tb. The resonant 
production of such a charged Higgs boson leads to a 
final state similar to the one resulting from single top 
s-channel production, qq — > W* — > tb. Therefore the 
analyses developed for the search for single top produc- 
tion can be applied to the search for a charged Higgs 
boson. Such an analysis was performed by the D0 collab- 
oration (|125[1 . in the topology arising from a subsequent 
t — > Wb — > Ivb decay. The large H+ mass, reconstructed 
from the decay products imposing the W boson and top 
quark mass constraints, was used as discriminating vari- 
able. No excess was observed over SM background pre- 
dictions, and upper limits were set on the production of 
a charged Higgs boson. The results are, however, not 
sensitive to Type II 2HDMs, but provide some exclusion 
in Type I 2HDMs. 
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V. SEARCHES FOR SUPERSYMMETRIC PARTICLES 
A. General features of SUSY models 

As explained in Sec. [TTJ, the main features of SUSY 
models for phenomenology are related to the type of me- 
diation mechanism for SUSY breaking, to the choice of 
soft breaking terms, and to whether or not i?-parity is 
assumed to be conserved. 

The most widely studied models involve gravity- 
mediation of SUSY breaking. In the minimal form of 
such models, mSUGRA, imparity is conserved, and only 
five parameters are needed beyond those already present 
in the standard model: a universal gaugino mass m\/2, a 
universal scalar mass mo, and a universal trilincar cou- 
pling A 0l all defined at the scale of grand unification, and 
tan/3 and the sign of /i. The low energy parameters, in- 
cluding are determined by the renormalization group 
equations and by the condition of electroweak symme- 
try breaking. In addition, it is commonly assumed that 
the LSP is the lightest neutralino Xi- A somewhat less 
constrained model keeps fi and vtia as independent low 
energy parameters, which is in effect equivalent to de- 
coupling the Higgs scalar masses from the masses of the 
other scalars. Such a model was largely used at LEP. 

Many studies have been performed where the assump- 
tion of i?-parity conservation is dropped, while keeping 
unchanged the other features of those mSUGRA inspired 
models. If i?-parity is violated, the superpotential is 
allowed to contain lcpton or baryon number-violating 
terms ff!26]) 

W Rp = \,,,l.,l.,lu, • A',. L,() ,1)-, ■ ,!),!),. , (39) 

where L and Q are lepton and quark doublet supcrficlds, 
E and D are lepton and down-type quark singlet super- 
fields, and i, j and k are generation indices. These terms 
are responsible for new couplings through which the LSP 
decays to SM particles. The simultaneous occurrence of 
different coupling types is however strongly constrained, 
e.g., by the bounds on the proton lifetime, which is why 
it is commonly assumed that only one of the i?-parity 
violating terms is present in the superpotential. 

In models with gauge-mediated SUSY breaking 
(GMSB), the LSP is a very light gravitino G, and the 
phenomenology is governed by the nature of the NLSP. 
In the minimal such model, mGMSB, all SUSY parti- 
cle masses derive from a universal scale A, and in most 
of the parameter space the NLSP is either the lightest 
neutralino \i or the lighter stau tr, the latter occur- 
ring preferentially at large tan (3. The couplings of the 
gravitino depend on yet another parameter, the SUSY- 
breaking scale y/F, which can be traded for the lifetime 
of the NLSP. 

Anomaly- mediation of SUSY breaking (AMSB) gener- 
ically leads to a neutralino LSP which is almost a pure 
wino W°, and has a small mass splitting with the lighter 
chargino. As a consequence, this chargino may acquire a 



phenomenologically relevant lifetime, possibly such that 
it behaves like a stable particle. 



B. Signatures and strategies 

Most of the searches for SUSY particles were per- 
formed within a "canonical scenario," the main fea- 
tures of which are borrowed from mSUGRA: i?-parity 
conservation, universal gaugino mass terms, a universal 
sfermion mass term, and a neutralino LSP. Because of 
i?-parity conservation, SUSY particles are produced in 
pairs, and each of the produced SUSY particles decays 
into SM particles accompanied by an LSP. Since the LSP 
is neutral and weakly interacting, it appears as missing 
energy, which is the celebrated signature of SUSY parti- 
cle production. 

Alternatively, if i?-parity is not conserved, the LSP 
decays to SM particles, so that no missing energy is ex- 
pected beyond that possibly arising from neutrinos. The 
signature of SUSY particle production is therefore to be 
sought in an anomalously large multiplicity of jets or lep- 
tons. The i?-parity violating couplings can also make it 
possible that SUSY particles are produced singly, rather 
than in pairs. 

In i?-parity conserving scenarios other than the canon- 
ical one, additional or different features are expected. In 
GMSB, each of the pair-produced SUSY particles decays 
into SM particles and an NLSP. The NLSP further decays 
into its SM partner and a gravitino. With a neutralino 
NLSP in the mass range explored up to now, the domi- 
nant decay is Xi ~^ jG, so that the final state contains 
photons, with missing energy due to the escaping grav- 
itinos. With a stau NLSP, the decay is tr — ► tG. If 
the stau lifetime is so long that it escapes the detector 
before decaying, the final state from stau pair production 
does not exhibit any missing energy, but rather appears 
as a pair of massive stable particles. A similar final state 
may also arise from chargino pair production in AMSB. 
Long-lived gluinos can lead to spectacular signatures if 
they are brought to rest by energy loss in the detector 
material. 

Except for the gluino, all SUSY particles are produced 
in a democratic way in e + e~ collisions via electroweak 
interactions. It is therefore natural that the searches at 
LEP were targeted toward the lightest ones. The re- 
sults of these searches could further be combined within a 
given model, thus providing constraints on the model pa- 
rameters. In contrast, it is expected that the most copi- 
ously produced SUSY particles in hadron collisions, such 
as pp at the Tevatron, will be colored particles, namely 
squarks and gluinos. Their detailed signature however 
depends on the mass pattern of the other SUSY particles, 
which may be present in their decay chains. This is why 
a specific model, usually mSUGRA, is needed to express 
the search results in terms of mass constraints. Thanks 
to lower masses and more manageable backgrounds, the 
search for gauginos produced via electroweak interactions 



17 



can be competitive at hadron colliders for model param- 
eter configurations where their leptonic decays are en- 
hanced. 

In e + e~ collisions, the production cross sections of 
SUSY particles are similar to those of their SM part- 
ners, except for the phase space reduction due to their 
larger masses. The data collected at the highest LEP 
energies, up to 209 GeV, are therefore the most relevant 
for SUSY particle searches. Mixing effects may however 
reduce these cross sections, as is the case for instance for 
neutralinos with a small Higgsino component, in which 
case the large integrated luminosity accumulated by the 
LEP experiments at lower energies also contributes to 
the search sensitivity. 

Although the center-of-mass energy of 1.96 TeV in 
pp collisions at the Tevatron allows higher new particle 
masses to be probed, large integrated luminosities are 
needed because of the rapid PDF fall off at high 
explained in Sec. |I.B| The search for SUSY particles at 
the Tevatron is also rendered more challenging than at 
LEP because of the large cross sections of the background 
processes. In the searches for squarks and gluinos, sig- 
nal production cross sections of the order of 0.1 pb at 
the edge of the sensitivity domain are to be compared to 
the total inelastic cross section of 80 mb. In the searches 
for gauginos, with similar signal production cross sec- 
tions in the mass range probed, the main backgrounds 
are W — > Iv and Z — > ££, with cross sections at the 
2.7 nb and 250 pb level per lepton flavor. 

In ep collisions at HERA, the most promising SUSY 
particle production process is single squark resonant pro- 
duction via an i?-parity violating X' 1 j 1 or A' llfe coupling, 
with a cross section depending not only on the squark 
mass, but also on the value of the coupling involved. The 
decay of the squark produced could be either direct, via 
the same A' coupling as for its production, or indirect 
through a cascade leading to the LSP, which in turn de- 
cays to two quarks and a neutrino or an electron. The 
mass reach at HERA is the full center-of-mass energy 
of 320 GeV, but the production of squarks with masses 
close to this bound involves quarks at large x values, so 
that the effective reach is substantially smaller, even for 
large values of the A' coupling. 



C. Searches in the canonical scenario 

As mentioned above, the characteristic signature of 
SUSY particle production in the canonical scenario is 
missing energy carried away from the detector by the 
LSPs at the end of the decay chains. 



1. Searches at LEP 

The main channels for SUSY particle searches in e + e~ 
coUisions are slepton (|T2"Tl HM), chargino (jT2^1 H"3tjj) 
and neutralino ([1311 11321 1133|) production. Squark pair 



production (|134p can also be relevant in some specific 
cases (|135p . 

Sleptons: In e + e~ annihilation, the search for SUSY 
particles that involves the least set of hypotheses for its 
interpretation is the search for smuons. Pair production 
proceeds via Z/j* exchange in the s-channel. Because 
of the small mass of the muon, the smuon mass eigen- 
states can be identified with the interaction eigenstates, 
of which jlft is the lighter one in models with slepton 
and gaugino mass unification. The search results were 
interpreted under this assumption, which is furthermore 
conservative, as the coupling of the fiR to the Z boson 
is smaller than that of the jl^. Only one parameter is 
needed to calculate the smuon pair production cross sec- 
tion, the smuon mass mp iR . The sole decay mode of a 
jjR NLSP is p,R — > fMXii so that smuon pair production 
leads to a final state consisting of two acoplanar muons 
with missing energy and momentum. The topology of 
this final state also depends on the mass of the LSP. If 
m^o is small, the final state is very similar to that arising 
from W pair production, with both W bosons decaying 
to a muon and a neutrino. If the pn — xl mass difference 
is small, the final state muons carry little momentum, so 
that the selection efficiency is reduced. In that configura- 
tion, the main background comes from "77 interactions," 
e + e~ — ► (e + )7*7*(e~) — * (e + )fi + p~ (e _ ), where the 
spectator electrons (e ± ) escape undetected in the beam 
pipe. The LSP mass to^o is therefore needed, in addition 
to the smuon mass, to interpret the search results. The 
constraints obtained in the (jrif tR , m^o ) plane by the four 
LEP experiments (136) are shown in Fig. 14 If the as- 



sumption that the smuon is the NLSP is dropped, further 
specification of the model is needed to turn the search 
results into mass constraints. An example is shown in 
Fig. [14] in the case of gaugino mass unification, for the 
specified values of p and tan (3. A slight reduction of the 
excluded domain is observed for low values of m,-,o , due 



to the competition of the fin — > pjQ. decay mode, with 
Depending on m^o, smuon masses smaller 



IXi ■ 



than 95 to 99 GeV are excluded, except for £ir — Xi mass 
differences below 5 GeV. 

Because of the larger r mass, compared to the muon 
mass, the hypothesis that the stau mass eigenstates can 
be identified with the interaction eigenstates may not 
hold, especially for large values of tan/3 that enhance 
the off-diagonal elements of the mass matrix in Eq. ( 29 1 . 



The coupling to the Z boson of the lighter stau mass 
eigenstate f\ may therefore be reduced with respect to 
the smuon coupling, and even vanish. Moreover, because 
there is at least one neutrino in each r decay, the visible 
energy of the final state arising from stau pair production 
is smaller than in the case of smuons, so that the selection 
efficiency is reduced. The mass lower limits obtained at 
LEP are therefore lower for staus than for smuons, from 
86 to 95 GeV, depending on m^o, provided the fi — xl 
mass difference is larger than 7 GeV (1136p . 

As for smuons, the selectron mass eigenstates can be 
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FIG. 14 Region in the (m^.m^o) plane excluded by the 
searches for smuons at LEP (136). The dotted contour is 
drawn under the assumption that the smuon decay branching 
ratio into fixi is 100%. 



identified with the interaction eigenstates. But because 
of the contribution of t-channel neutralino exchange to se- 
lectron pair production, the gaugino sector of the model, 
mass spectrum and field contents, has to be specified to 
interpret the results of the searches for acoplanar elec- 
trons. With gaugino mass unification and for tan j3 = 1.5 
and \x — —200 GeV, a selectron mass lower limit of 
100 GeV was obtained for m^o < 85 GeV <fiM\i . Neu- 
tralino t-channcl exchange can furthermore mediate as- 
sociated e^efl production. This process is useful if the 
e.R — Xi mass difference is small, because the electron 
from the — > ex\ decay can be energetic enough to lead 
to an apparent single electron final state. Both gaug- 
ino and slepton mass unifications have to be assumed 
for the masses of the two selectron species to be related. 
Under these assumptions, a lower limit of 73 GeV was 



me R , independent of the 



X? 



mass differ- 



set on 
ence 

From the measurement of the invisible width of the Z 
boson (flU)l , a general mass lower limit of 45 GeV can be 
deduced for a sneutrino LSP or NLSP. 

Charginos and neutralinos: As evident from 
Eq. ( 28 1 , three parameters are sufficient to fully spec- 



ify the masses and field contents in the chargino sector. 
These may be taken to be M 2 , /i, and tan (3. The lighter 
of the two charginos will simply be denoted "chargino" 
in the following. To specify the neutralino mass ma- 
trix of Eq. (26 1, one more parameter, Mi, is needed. 



If gaugino mass unification is assumed, the two gaugino 
masses are related by Mi = (5/3) tan 2 Qy/Mi ~ O.5M2. 



Unless otherwise specified, this relation is assumed to 
hold in the following. Charginos are pair produced via 
s-channcl Z/7* and t-channel v e exchanges, the two pro- 
cesses interfering destructively. The three-body final 
states ff'xi are reached in chargino decays via virtual 
W or sfermion exchange. If kinematically allowed, two- 
body decays such as x ~ > are dominant. Simi- 
larly, neutralino pair or associated production proceed 
via s-channel Z and t-channel selectron exchanges, and 
X2 three-body decays to ffxi y i a virtual Z or sfermion 
exchange; whenever kinematically allowed, two-body de- 
cays such as x§ — * y v are dominant. 

If sfermions are heavy, chargino decays are mediated by 
virtual W exchange, so that the final states arising from 
chargino pair production are the same as for W pairs, 
with additional missing energy from the two neutralino 
LSPs: all hadronic (qq' XiQQ' Xi) , mixed (qq'xi^Xi), an d 
fully leptonic (i^Xi^Xi)- Selections were designed for 
these three topologies and for various m^± — m^o regimes, 
with no excess observed over SM backgrounds. From a 
scan over M2, At, and tan/3, a chargino mass lower limit 
of 103 GeV was derived for m v > 200 GeV ££39). For 
smaller sneutrino masses, the limit is reduced by the de- 
structive interference in the production. This limit holds 
for Mi <~ 1 TeV. For larger M2 values, the selection 
efficiency decreases rapidly as the x ± — X? mass differ- 
ence becomes smaller. If this mass difference becomes so 
small that even the x ~* 7r± Xi decay mode is closed, the 
chargino becomes long lived. Searches for charged mas- 
sive stable particles, in which advantage is taken of their 
larger ionization power, were designed to cope with this 
configuration. For slightly larger mass differences, the 
visible final state is so soft that even triggering becomes 
problematic. Chargino pair production can however be 
"tagged" by an energetic photon from initial state radi- 
ation, e + e~ — > 7X + x~, providing access to those almost 
invisible charginos, although at a reduced effective cen- 
ter of mass energy. The combination of these analysis 
techniques allowed chargino masses smaller than 92 GeV 
to be excluded, irrespective of the x^ 1 ~ Xi mass differ- 
ence dnoj. 

For lower sfermion masses, the sensitivity of the for- 
mer analyses is reduced first because of the destructive 
interference between the s-channel Z/j* and t-channel 
sneutrino exchanges, and second because of the open- 
ing of two-body decays. The latter effect is specifically 
detrimental in the "corridor" of small — v mass dif- 
ferences, where the final state from the x — * tv de- 
cays becomes invisible in practice. Gaugino mass unifi- 
cation allows indirect limits on charginos to be obtained, 
based on constraints on the parameter space resulting 
from searches for pair or associated neutralino produc- 
tion, e.g., e + e~ — > X2X2 or XiX°- I n order to relate 
all production cross sections and decay branching frac- 
tions, it is however necessary to fully specify the sfermion 
spectrum, which is done with the assumption of sfermion 
mass unification. The results of the chargino and neu- 
tralino searches are then expressed as exclusion domains 
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FIG. 15 Regions in the (/it, M%) plane excluded by the LEP I 
constraints (1), and by the searches for charginos (2), neu- 
tralinos (3) and sleptons (4) at LEP II, for tan/3 = V2 and 
mo = 75 GeV. The region (5) is excluded by the Higgs boson 
searches at LEP II. This figure is from Ref. I|14ip . 
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FIG. 16 Lower mass limit for the lightest neutralino as a 
function of tan/3, inferred in the conventional scenario from 
searches at LEP for charginos, sleptons, and neutral Higgs 
bosons (|142l) . The dashed contour is the limit obtained for 
large mo- 



rn the (/i, M-i) plane for selected values of tan/3 and 



mo. The invisible two-body decay X2 



can how- 



ever cause a large sensitivity reduction. Since this con- 
figuration occurs for low Too values, constraints arising 
from the slepton searches can be used to mitigate this 
effect. With gaugino and sfermion mass unification, the 
slepton masses are related to the model parameters by 
to?- ~ m 2 , + 0.22Aff — sin 2 9 w m 2 z cos 2/3, so that a limit 
on rn^ R can be turned into a limit on M2 for given val- 
ues of tan (3 and too- After a proper combination of the 
searches for charginos, neutralinos and sleptons, an ex- 
ample of which is shown in Fig. |15[ it turns out that 
the chargino mass limit obtained in the case of heavy 
sfermions is only moderately degraded. 

Direct searches for the lightest neutralino had been 
performed at lower energy e + e _ colliders, PEP and PE- 
TRA, in the reaction e + e~ — > 7X1X11 where the photon 
from initial state radiation is used to tag the production 
of an invisible final state. At LEP, at or above the Z 
resonance, the irreducible background from e + e _ — > 
is too large to obtain competitive results. Furthermore, 
production via s-channel Z exchange may simply vanish, 
e.g., if the LSP is a pure photino, while production via 
^-channel selectron exchange can be made negligible if se- 
lectrons are sufficiently heavy. Indirect limits on the mass 
of the LSP can however be obtained within constrained 
models. With gaugino mass unification, to^o is typically 
half the chargino mass. As a result, the chargino mass 
limit translates into a x? mass lower limit of 52 GeV for 



heavy sfermions and large tan [3. If sfermion mass unifi- 
cation is used in addition, a limit of 47 GeV is obtained 
at large tan/3, independent of mo. This limit is set by 
searches for sleptons in the corridor. For low values of 
tan/3, constraints from the Higgs boson searches can be 
used, as was shown in Sec. |III.B| for benchmark scenar- 
ios. A complete scan over too, toi/ 2 , M an d tan/3 was 
performed and, for each parameter set, the maximum h 
mass predicted was compared to the experimental limit, 
and the constraints from chargino and slepton searches 
were included. The translation of the scan result in terms 
of excluded domain in the (tan /3, to-^o) plane is shown in 
Fig. 16 from which a Xi mass lower limit of 47 GeV is 



derived (|142|l . Within the more constrained mSUGRA 
scenario, wherein /1 is calculated from the other parame- 
ters, this limit becomes 50 GeV (|143p . 

Squarks: On general grounds, the mass reach for 
strongly interacting particles is expected to be substan- 
tially higher at the Tevatron than at LEP. For some spe- 
cific configurations, however, the searches at the Teva- 
tron become inefficient, in which cases the results ob- 
tained at LEP remain of interest. This is particularly 
relevant for third generation squarks which may be sub- 
stantially lighter than the other squarks, as motivated 
in Sec. [TTJ The lighter third generation mass eigenstates 
are simply denoted stop and sbottom, t and b, in the 
following. 

In the mass range accessible at LEP, and given the 
chargino mass limit which effectively forbids t — > b\ + , the 
stop is expected to decay into a charm quark and a neu- 
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tralino, t — > cxf, as long as mj < miy + mj + m^o (|144p . 
Because this decay is a flavor-changing loop process, the 
stop lifetime can be large enough to compete with the 
hadronization time, and the simulation programs were 
adjusted to take this feature into account. The final 
state from stop pair production exhibits an acoplanar jet 
topology, for which no signal was observed above stan- 
dard model backgrounds. As already explained for staus, 
the amount of mixing between the weak cigcnstatcs can 
be such that the stop does not couple to the Z boson. 
In this worst case scenario, stop mass lower limits rang- 
ing from 96 to 99 GeV were obtained, depending on the 
Xi mass, as long as — m^o — m c > 5 GeV (|145p . For 
smaller t — xi mass differences, long-lived i?-hadrons may 
be produced in the stop hadronization process. The pro- 
duction of such i?-hadrons and their interaction in the 
detector material were taken into account in a dedicated 
search, from which a stop mass lower limit of 63 GeV was 
derived, valid for any rrij: — to^o (|146[l . For specific pa- 
rameter choices, and in spite of the slepton mass limits, 
it can be that the t — > b£D decay is kinematically allowed, 
in which case it is dominant. From a search for events 
exhibiting jets, leptons and missing energy, a stop mass 
lower limit of 96 GeV was obtained, valid for sneutrino 
masses smaller than 86 GeV (|145p . 

The case of a light sbottom is much simpler, as the tree- 
level b — > b\\ decay mode is dominant. From searches 
for acoplanar 6-flavored jets, a mass lower limit of about 
95 GeV was obtained in the worst case scenario where 
the sbottom does not couple to the Z (j!45p . 



2. Searches at the Tevatron 

The program most widely used for the calculation of 
SUSY particle production cross sections at the Tevatron 
is PROSPINO ()147p . which provides next-to-leading order 
accuracy. The results reported below were generally ob- 
tained with the CTEQ6.1M PDF set (Hi [MD- Var- 
ious codes were used to calculate the low energy SUSY 
spectrum from initial parameters at the grand unification 
scale: SUSPECT flM]), SOFTSUSY (52), and isajet (|5T]l . 
This may introduce slight inconsistencies when compar- 
ing results in different channels or from different experi- 
ments in terms of parameters at the high scale. The pro- 
duction of SUSY particles was in general simulated with 
PYTHIA (|150p . with decays modeled with SDECAY (|15ip 
or with ISASUGRA as implemented in pythia. Typically, 
SM backgrounds were simulated with ALPGEN (I152p for 
the production of W and Z bosons in association with 
jets, or with PYTHIA otherwise. 

As already mentioned in Sec. |V.B[ the channels most 
relevant for SUSY particle searches at hadron colliders 
are the production of squarks and gluinos on the one 
hand, of electroweak gauginos on the other. For squarks 
and gluino, the search is conducted in events exhibiting 
a jets+|?T topology (1531 11541 1155p . while for electroweak 
gauginos, it is conducted in the trilepton final state (|156l 



HSU USB). 

Generic squarks and gluinos: Depending on the 
squark and gluino mass hierarchy, different pair produc- 
tion processes via the strong interaction are expected to 
contribute in pp collisions at the Tevatron: qq and, to 
a lesser extent, qq, if rriq <C vn~ g \ gg if nig <C rriq; and 
all of these processes, as well as qg, if the squark and 
gluino masses are similar. If rriq < m~ g , squarks are ex- 
pected to decay directly into a quark and a gaugino, 
q — > qx, where x is dominantly xi for q^, and x or 
X® for Ql- If m g < m qi gluinos are expected to decay 
via virtual squark exchange into a quark, an antiquark, 
and a gaugino, g — > qqx, where x 1S typically or 
X° l|159j I160p . The heavier gauginos further decay into a 
fermion-antifcrmion pair and an LSP, Xi i so that there is 
always some missing Et in the final state. More detailed 
predictions can be made only within a specific model such 
as mSUGRA. 

The aforementioned production processes have been 
searched for by CDF and D0 in topologies involving at 
least two jets, four jets and three jets, all with larger- 
Initial and final state radiation of soft jets can increase 
further those jet multiplicities. A first class of back- 
ground to squark and gluino production arises from pro- 
cesses with intrinsic , such as (W — > £v)+iets, where 
the lepton escapes detection, or (Z — > w)+jets, which 
is irreducible. Monte Carlo simulations were used to 
estimate those backgrounds, after calibration on events 
where leptons from W — > Iv or Z — > 11 are detected. An- 
other class of background is due to multijet production 
by strong interaction. Although there is no intrinsic IpT 
in such events, fake^^ can arise from jet energy mismea- 
surements (and also real^r from semileptonic decays of 
heavy flavor hadrons). In such events, the^V distribu- 
tion decreases quasi-exponentially, and the direction of 
the JfjT tends to be close to that of a mismeasured jet. 
Requiring sufficiently large JpT and applying topological 
selection criteria allows this background to remain un- 
der control. While D0 applied criteria tight enough to 
reduce this background to a negligible level, CDF esti- 
mated its remaining contribution based on simulations 
calibrated on control samples. 

No excesses of events were observed over SM back- 
grounds, which was translated into exclusion domains in 
the plane of squark and gluino masses. To this end, a spe- 
cific SUSY model had to be chosen, so that the masses 
and decay modes of all the gauginos involved in the de- 
cay chains could be determined. The model used by both 
CDF and D0 was mSUGRA, with A = 0, fx < 0, and 
tan P = 5 (CDF) or 3 (D0). The production of all squark 
species was considered, except for the third generation 
(CDF) or for stops (D0), and the squark mass quoted 
was the average of the masses of the squarks considered. 
Finally, the large theoretical uncertainties associated to 
the choices of PDFs and of the factorization and renor- 
malization scales had to be taken into account when turn- 
ing cross section upper limits into exclusion domains in 
terms of masses. Based on an integrated luminosity of 
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FIG. 17 Region in the (m,g,mq) plane excluded by D0 (j!61|) 
and by earlier experiments. The red curve corresponds to the 
nominal scale and PDF choices. The yellow band represents 
the uncertainty associated with these choices. The blue curves 
represent the indirect limits inferred from the LEP chargino 
and slepton searches. 



2.1 fb _1 , D0 excluded the domain shown in Fig. 17 from 
which lower limits of 379 and 308 GeV were derived for 
the squark and gluino masses, respectively, as well as a 
lower limit of 390 GeV if rriq = rrig (|161| l. Similar results 
were obtained by the CDF collaboration ([162ft . 

Third generation squarks: As already mentioned, 
a stop NLSP decays into a charm quark and a neutralino 
as long as < mw+mb+m^o . The final state from stop 
pair production therefore consists in acoplanar charm jets 
and Tfij- . Because only one of the squark species is now 
produced, the cross section is smaller than for generic 
squarks, and the mass reach is therefore lower. As a 
consequence, the jets are softer, and there is also less 
IpT ■ The corresponding loss of sensitivity was attenuated 
by making use of heavy-flavor tagging, which resulted 

obtained by 
of data. It can 



in the exclusion domain shown in Fig. 
D0 (|163p from an analysis of 1 fb _1 
be seen that a stop mass of 150 GeV is excluded for 
to^o = 65 GeV. In spite of the larger mass reach at the 
Tevatron, the LEP results remain the most constraining 
for t— Xi mass differences smaller than ~ 40 GeV. Similar 
searches were performed for a sbottom NLSP decaying 
into b\i (|164l 1165ft , with better sensitivity due to a more 
efficient heavy- flavor tagging for b than for c quarks. A 
mass lower limit of 222 GeV was obtained by D0 for 
m^a < 60 GeV, based on 310 pb -1 of data. 

Other mass hierarchies were considered, where the stop 
or sbottom is not the NLSP. Three-body stop decays, 
t — > b£v, are dominant if kinematically allowed and when 
t — > x + b is not, which is possible for some model param- 
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FIG. 18 Region in the (mj, m^o) plane excluded by D0 i 
and by earlier experiments. The solid curve corresponds to 
the nominal scale and PDF choices. The yellow band repre- 
sents the uncertainty associated with these choices. 



eter choices in spite of the mass limits on charged slep- 
tons available from LEP. The final states investigated by 
D0 comprised two muons or a muon and an electron, 
with b jets and-^r- Based on an analysis of 400 pb -1 
of data, the largest stop mass excluded is 186 GeV, for 
rn,y = 71 GeV (166). If the chargino is lighter than the 
stop, the t — > bx + decay is dominant. A search was per- 
formed by CDF in the two lepton, two b jets and^r final 
state, with a sensitivity depending on the branching frac- 
tion of the chargino leptonic decay, x^ ~ * ^Xij which is 
enhanced for light sleptons. An example of an excluded 
domain in the (m t -, m^o) plane is shown in Fig. 19 (167), 
based on 2.7 fb _1 of data. In both of those searches, 
the background from top quark pair production was a 
major challenge. Yet another mass hierarchy was consid- 
ered by CDF, namely that where the sbottom is the only 
squark lighter than the gluino. In such a configuration, 
the g — > bb decay is dominant, and gluino pair produc- 
tion then leads to a final state of four b jets and^r- This 
search was performed in a data sample of 2.5 fb _1 , and 
lead to excluded sbottom masses as large as 325 GeV 
for gluino and LSP masses of 340 and 60 GeV, respec- 
tively (fT68jl . 

Charginos and neutralinos: The associated pro- 
duction of charginos and neutralinos, pp — > X^xH, is 
an electroweak process mediated by s-channel W and 



t-channel squark exchanges. Leptonic decays, x^ ~ * 
£ i/Xi an d X2 ~ ¥ £ + £~Xii are mediated by W and Z ex- 
change, respectively, and by slepton exchange. If sleptons 
are light, leptonic decays can be sufficiently enhanced for 
searches in final states consisting of three leptons and^r 
to become sensitive in spite of production cross sections 
of a fraction of a picobarn. An additional challenge is the 
rather small energy carried by the final state leptons in 
the chargino and neutralino mass domain to which the 
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FIG. 19 Regions in the (mj, m-o) plane excluded by 
CDF CUT} for m i± = 125.8 GeV and for various values of 
the branching fraction for the x — » ifXi decay. 
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FIG. 20 Regions in the (mo,mi/2) plane excluded by the 
D0 search for trileptons l)170p . 



D. Searches in non-canonical scenarios 



searches at the Tevatron are currently sensitive. 

In both the CDF ff!69]) and D0 (|T70j) analyses, only 
two leptons were required to be positively identified as 
electrons or muons. 4 Allowing the third lepton to be 
detected as an isolated charged particle track provided 
sensitivity to final states including a r lepton that de- 
cays into hadrons. In the CDF analysis, the trilepton fi- 
nal state was split into topologies with different signal to 
background ratios, the purest being when the three lep- 
tons are positively identified as electrons or muons with 
tight criteria. In the D0 approach, different selections 
were optimized according to the amount of energy avail- 
able to the lepton candidates. The ultimate background 
for these trilepton searches is associated WZ production. 

The D0 search, based on an integrated luminosity of 
2.3 fb -1 , excludes regions in the mSUGRA parameter 
space as shown in Fig. [20] for Aq = 0, tan/3 = 3 and 
/i > 0. It can be seen that the domain excluded at LEP 
is substantially extended by these trilepton searches. The 
interruption in the exclusion domain is due to configura- 
tions where the small x® — ^ mass difference results in one 
of the final state leptons carrying too little energy, thus 
preventing efficient detection. Requiring only two lep- 
tons to be identified, but with same charge sign in order 
to reduce the otherwise overwhelming SM backgrounds, 
should provide sensitivity in that region, as was shown by 
D0 in an analysis based on a smaller data sample (171). 
The same-sign dilcpton signature was also considered in 
an earlier CDF analysis ()172p . 



4 The D0 analysis also considered final states with a muon and 
one or two r leptons identified. 



1. i?-parity violation 

Searches for SUSY with i?-parity violation were per- 
formed at LEP, the Tevatron and HERA. Both i?-parity 
conserving pair production of SUSY particles and R- 
parity violating resonant single SUSY particle produc- 
tion were considered. The produced particles were subse- 
quently subject to either direct or indirect (via a cascade 
to the LSP) i?-parity violating decays. Unless otherwise 
specified, a single i?-parity violating coupling is assumed 
to be non- vanishing in the following, large enough for the 
lifetime of the LSP to be safely assumed to be negligible. 

Searches at LEP: Extensive searches for pair pro- 
duction were performed at LEP, involving all possible 
i?-parity violating couplings. The possible final states 
are numerous, ranging from four leptons and missing en- 
ergy for Xi P a i r production, with decays mediated by a 
A-type coupling, e.g., Xi ~^ e^iu, to ten hadronic jets 
and no missing energy for chargino pair production, with 
X — > qq'xi followed by a x? decay into three quarks via 
a A"-type coupling, e.g., Xi — > u dd. The results of these 
searches are at least as constraining as in the canonical 
scenario (11751 [T71 H751 ITTrJl IT771) . 

The production of a sneutrino resonance via a Ayi 
coupling was also investigated. No signal was observed, 
and mass lower limits almost up to the center-of-mass en- 
ergy were set for sufficiently large values of the i?-parity 
violating coupling involved (|T751 ffffl USUI H5T1IT52]! . 

Searches at HERA: As explained in Sec. |V.B[ 
the HERA ep collider is most effective in the searches 
for i?-parity violating resonant single squark produc- 
tion via a A'-type coupling. Direct and indirect squark 
decays were investigated, and the search results were 
combined to lead to squark mass lower limits up to 
275 GeV (IT831 HMj ITSS) . within mild model assump- 
tions, for a A' coupling of 0.3, i.e., with electromagnetic 
strength. 
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Searches at the Tevatron: A fully general search 
for all i?-parity violating couplings is not possible at the 
Tevatron, as it was at LEP. For instance, A" couplings 
lead to multijct final states with no or little missing en- 
ergy, which cannot be distinguished from standard mul- 
tijet production. Searches have therefore been designed 
for specific choices of couplings leading to distinct signa- 
tures. 

Gaugino pair production followed by indirect decays 
has been extensively studied by both CDF f|186|) and 
D0 (|187p in the case of a A-type coupling. The final state 
is expected to contain four charged leptons, with flavors 
depending on the indices in the Xijk coupling, and^r due 
to two neutrinos. For mo = 1 TeV, tan /3 = 5, and /i > 0, 
the chargino mass lower limits obtained by D0 from an 
analysis of 360 pb _1 of data are 231, 229, and 166 GeV 
for the A121, A122, and A133 couplings, respectively, with 
reduced sensitivity in the last case due to the occurrence 
of t leptons in the final state. 

Stop pair production, with t — ► br via a A 333 coupling 
has been searched by CDF (|188p in the topology where 
one r lepton decays into an electron or a muon, and the 
other into hadrons. From an analysis of 322 pb _1 of data, 
a stop mass lower limit of 151 GeV was derived. 

Resonant smuon or sneutrino production could be me- 
diated by a A 211 coupling. With indirect decays, the 
final state would exhibit at least one muon and two jets. 
This topology was investigated by D0 (|189p . and an ex- 
cluded domain was set in the (m^,A 211 ) plane, leading 
to a smuon mass lower limit of 363 GeV for A 211 = 0.1, 
and for A n = 0, tan j3 = 5, and // < 0. 

Resonant sneutrino production mediated by a A^ 1:L cou- 
pling was also investigated by CDF and D0 (fl"901 [T9T1 
11921 I193p . now assuming that the sneutrino decays di- 
rectly via a A-type coupling. The final states considered 
were ee, e/z, /z/i, and rr. The sneutrino mass limits ob- 
tained depend on the product of the two couplings in- 
volved. 



2. Gauge-mediated SUSY breaking 

As already explained, the LSP in GMSB is a very light 
gravitino, and the phenomenology depends essentially on 
the nature of the NLSP, a neutralino or a stau, possibly 
almost mass degenerate with and jlR, and on its life- 
time. 

Neutralino NLSP: In the mass range of current in- 
terest, a neutralino NLSP decays into a photon and a 
gravitino, xl ~ * lG. Pair production of such a neutralino 
at LEP would therefore lead, assuming prompt decays, 
to a final state of two acoplanar photons and missing en- 
As can be seen in Fig. 
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no excess was observed 
above the SM background from e + e _ — > (Z^*) — > vv)^. 
In GMSB, Xi nas a large Bino component, so that 
pair production in e + e~ interactions proceeds via se- 
lectron i-channcl exchange. An excluded domain in the 
(ms R , rn^o) plane was therefore derived (] 194f) . ruling out 
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FIG. 21 Mass of the invisible system recoiling against pairs 
of photons at LEP (TngJ). 



the GMSB interpretation (in terms of selectron pair pro- 
duction) of an anomalous eejj+ fx event that had been 
observed by CDF (195 ) during Run I of the Tevatron. 

Searches were also performed at LEP for photons not 
pointing toward the interaction point, which could arise 
from non-prompt decays of a neutralino NLSP. For even 
longer lifetimes, the phenomenology becomes identical 
to that of the canonical scenario. The results of these 
various searches for a neutralino NLSP were combined 
with those in various topologies expected to arise from 
heavier SUSY particle production to lead to a robust 
neutralino mass lower limit of 54 GeV within the minimal 
GMSB framework ([155111571) . 

Searches for acoplanar photons with large J^t were 
performed at the Tevatron by both CDF (198) and 
D0 (p.99). This topology is expected to arise when- 
ever SUSY particles are pair produced, which subse- 
quently decay to a neutralino NLSP with negligible life- 
time. No excess of events was observed over the back- 
grounds due to photon misidcntification or from fake^r, 
all determined from data. These results were interpreted 
within the "Snowmass slope SPS 8" benchmark GMSB 
model (200J where the only free parameter is the effective 
SUSY breaking scale A. The other parameters were fixed 
as follows: 7V 5 = 1 messenger, a messenger mass of 2 A, 
tan/3 = 15, and fi > 0. Neutralino NLSP masses smaller 
than 138 GeV are excluded by the CDF analysis, based 
on 2 fb _1 of data. 

The possibility of non-prompt neutralino NLSP decays 
was also investigated by CDF ()20ip . making use of the 
timing information of their calorimeter. No signal of de- 
layed photons was observed in a data sample of 570 pb _1 , 
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NLSP mass and lifetime (1951 l20lj) . 



from which an excluded domain in the plane of the mass 
and lifetime of the NLSP was inferred, as shown in Fig. [22] 
together with the result of Ref . (|198[l . 

Stau NLSP: For prompt f — > tG decays, the fi- 
nal state arising from stau pair production at LEP is 
the same as in the canonical scenario with a very light 
Xi- For very long lifetimes, the searches for long lived 
charginos already reported apply. Searches for in-flight 
decays along charged particle tracks were designed to ad- 
dress intermediate lifetimes. The combination of all these 
searches allowed a stau NLSP mass lower limit to be set 
from 87 to 97 GeV, depending on the stau lifetime, as 
shown in Fig. [23] (|2"02"]). 



3. Other non-canonical scenarios 

A number of searches were performed at LEP and at 
the Tevatron in other non-canonical scenarios. 

Stable charged particles: In anomaly mediated 
SUSY breaking, the LSP is wino-like, and the x — Xi 
mass difference is therefore small. As a result, stable 
charginos are not unlikely. The searches designed at LEP 
in the canonical scenario for large Mi values apply here 
equally well. At the Tevatron, a search was performed 
by D0 for pairs of charged massive stable particles that 
could result from chargino pair production. Such parti- 
cles would behave like slow moving muons that could be 
detected as delayed signals in the muon system. No sig- 
nificant excess of delayed muons was observed in 1.1 fb -1 
of data, and a mass lower limit of 206 GeV was set on 
long-lived wino-like charginos (203). 

A search for stable stops was performed by CDF in 
1 fb -1 of data, using a high pr muon trigger and their 
time-of-flight detector. Stable stops hadronize to form 
i?-hadrons which behave like slow muons. A model for 




FIG. 23 Domain excluded at LEP (f2"02|) in the (log 10 r, m f ) 
plane, where r is the lifetime in seconds and rrif is the mass 
of a stau NLSP in GMSB. The green shaded region is the 
excluded domain, and the dashed blue curve is the expected 
boundary of the exclusion region. 



the interactions of those i?-hadrons with the detector ma- 
terial was constructed, within which a stop mass lower 
limit of 249 GeV was derived flHI|). 

Stable or long-lived gluinos: Models have been 
built where the gluino could be the LSP and therefore 
stable, if imparity is conserved (|2U5l I2UB1 |2U7| . Al- 
ternatively, gluinos may decay, but with long lifetimes. 
This occurs, for example, in models with "split SUSY," 
unnatural models in which all squarks and sleptons are 
very heavy, but the gauginos remain at the electroweak 
scale (208; 209i I210p . Since gluino decays are mediated 
by squark exchange, the gluino becomes long-lived. 

Although gluinos cannot be produced directly in e + e~ 
interactions, they could be produced via gluon split- 
ting, e.g., e + e~ — > qqg* — > qqgg, and hadronize into 
metastable "i?-hadrons." The QCD predictions for four- 
jet events would therefore be modified. Gluinos could 
also be produced in the decay of heavier squarks. Dedi- 
cated analyses were performed at LEP (|211II212"]| . leading 
to a mass lower limit of 27 GeV for a stable gluino. 

A search for long-lived gluinos was also performed by 
D0 with 410 pb^ 1 of data (|213p . After hadronization 
into an i?-hadron, a long-lived gluino could come to rest 
in the calorimeter and decay later on, during a bunch 
crossing different from the one during which it was cre- 
ated (|214[1 . The main decay mode expected is g — > gx\, 
leading to an hadronic shower originating from within the 
calorimeter and not pointing toward the pp interaction 
region, in an otherwise empty event. No excess of this 
anomalous topology was observed over the background 
due to cosmic muons or to the beam halo. The gluino 
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mass lower limits derived depend on the lifetime Tg, the 
branching fraction B for the decay mode considered, the 
Xi mass, and the cross section or for the conversion of a 
neutral i?-hadron into a charged one in the calorimeter. 
As an example, a mass limit of 270 GeV was obtained 
for Tg < 3 hours, B = 1, rn^a = 50 GeV, and ctr = 3 mb. 



VI. SUMMARY 

Supersymmetry is one of the most promising ideas for 
extending the standard model. When realized at the 
weak scale, many SUSY models provide natural, and 
even elegant, solutions to the most pressing problems 
in particle physics today by stabilizing the gauge hier- 
archy, providing dark matter candidates, and accommo- 
dating force unification, both with and without gravity. 
In addition, the general framework of weak-scale SUSY 
is flexible enough to encompass a wide variety of new 
phenomena, including extended Higgs sectors, missing 
energy, long-lived and metastable particles, and a host 
of other signatures of new physics. Searches for SUSY 
are therefore also searches for other forms of new physics 
which, even if less profoundly motivated, are, of course, 
also important to pursue. 

In this review, we have comprehensively summarized 
the state of the art in searches for SUSY at the high 
energy frontier. Although this is a continuously evolving 
subject, this review provides a snapshot of the field at a 
particularly important time, when final results from LEP 
and HERA arc in hand, the Tevatron experiments have 
reported deep probes of many supersymmetric models 
with several fb _1 of data, and the LHC will soon begin 
operation. 

This review has summarized searches for both super- 
symmetric Higgs bosons and standard model superpart- 
ners. In the Higgs sector, SUSY requires a light neutral 
Higgs boson. This Higgs boson could be standard-model 
like, but it could also have non-standard couplings. In 
addition, it is accompanied by other Higgs bosons, both 
neutral or charged. The most stringent constraints on a 
SM-like Higgs boson currently come from LEP, with a 
mass lower limit of 114.4 GeV that applies in the MSSM 
at low tan/3. Furthermore, the LEP experiments set a 
lower limit of 93 GeV on the lightest neutral Higgs boson 
of the MSSM, independent of tan 0. The MSSM parame- 
ter space has now been further restricted by the Tevatron 
experiments. For example, tan j3 values larger than 40 are 
excluded for = 140 GeV. For charged Higgs bosons, 
LEP excludes masses below 78.6 GeV, and the Tevatron 
experiments have extended this mass limit to ~ 150 GeV 
for very large values of tan (3. 

For superpartners, the bounds are, of course, model- 
dependent, but the main results may be summarized 
as follows. The searches at LEP have constrained the 
masses of all SUSY particles, except for the gluino and 
the LSP, to be larger than approximately 100 GeV in 
most SUSY scenarios. Furthermore, an indirect lower 



limit on the mass of a neutralino LSP has been set 
at 47 GeV in the MSSM with gaugino and sfermion 
mass unification. The higher center-of-mass energy at 
the Tevatron has allowed tighter mass limits to be ob- 
tained for strongly interacting SUSY particles: 379 and 
308 GeV for squarks and gluinos, respectively, within 
the mSUGRA framework at low tan (3. In that same 
model, domains beyond the LEP reach were also probed 
by searches for associated chargino-neutralino produc- 
tion. 

In the near future, the first indication for SUSY at 
high energy colliders could be the observation of a light 
neutral Higgs boson at the Tevatron. Of course, such 
a discovery is not proof of SUSY — only the discovery 
of superpartners would unambiguously establish SUSY 
as being realized in nature. Once collisions begin at the 
LHC and the detectors are sufficiently understood, it will 
not take more than ~ 1 fb _1 to discover squarks and 
gluinos with masses less than ~ 1.5 TeV ([2T1 |2"3")) . A 
new era will then begin during which the whole SUSY 
spectrum will have to be deciphered, and the properties 
of the SUSY model established. Many more fb _1 will be 
needed for that purpose, and to unravel the spectrum of 
SUSY Higgs bosons. 
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